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NOTATION 


Ihe nomenclature used in this report agrees generally with that recommended 
in "Soil Mechanics Nomenclature” , ASCS Manual of Engineering Practice No. 22, 
19Ul. A partial list of recommended symbols and new symbols is given below: 


Symbol Description 

a ; , Thermal Diffusivity 

C s Specific Heat of Soil Solids 

Cf ■“ Volumetric Heat Capacity of Frozen Soil 

Cu Volumetric Heat Cauacitv of Unfrozen 

Soil 

F Freezing Index 

kf Thermal Conductivity of Frozen Soil 

k u Thermal Conductivity of Unfrozen Soil 

L Latent Heat of Fusion of Soil Moisture 

n Surface Transfer Coefficient 

q Rate of Heat Flow 

As Finite Distance 

t Time 

(For frost penetration formulas and 
nomograph: Duration of Freezing Index) 

u Internal Energy, total heat content 

v Temperature 

v Degrees F by which mean annual temper - 

^ ature exceeds freezing point of soil 

moisture 

w 'Water Content of Soil 

iii 


Units 
sq.ft, per hr. 

BTU per lb. 

BTU per (cu.ft.X°F) 

BTU per (cu.ft.)( °F) 
Degree Days 
BTU per (ft. )(°F)(hr.) 
BTU per (ft.)( °F)(h r.) 
BTU per cu.ft. 
(Dimensi onless ) 

BTU per (sq.ft .)( hr. ) 

ft. 

hr. 

(days ) 

3TU per cu.ft. 

°F 

°F 
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Symbol 


Units 


i . . 

Description 


X 

Ax,4y,A2 


(X 


Depth of Frost Penetration 

Finite Distances in the x, y, and 
z directions 

Thermal Ratio 

a A t 

- As 2 

Correction Coefficient 
Fusion Parameter 


feet (or inches) 
ft. 

(Dimensionless ) 

(Dimensionless ) 

(Dimensionless ) 
(Dimens i onles s ) 
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SYNOPSIS 

I * 

This report contains the results of studies related to the prediction 
of the frost penetration depth below ground surfaces* In addition, the so- 
called thaw-consolidation problem has been studied and the results are repor- 
ted by the writers.- 

Part II is devoted to the presentation of a mathematical formulation of 
the problem derived from first principles of heat transfer. The fluid and 
electric analogies to thermal problems are treated in some detail since they 
are of especial value not only for visualizing the thermal problem but as 
important aids in solving complex situations. Furthermore, because a closed- 
form analytical solution is impossible except for the simplest cases, consi- 
derable effort has been devoted to a discussion of practical methods of com- 

■a 

putations such as numerical and machine solutions. Descriptions of an hy- 
draulic model and an electronic analog computer applicable to freezing and 
thawing problems are presented. 

The second section of Part II deals with the writers' rational formula 
for the depth of frost penetration. Results of statistical studies involving 
actual frost penetration data are described. The effects of cyclic variations 
in surface temperature and of multilayered systems on the depth of freezing 
are presented. 

Part III deals with all phases of the writers' depth of frost penetration 
nomograph including instructions for its use and its limitations and assump- 
tions. A copy of the nomograph is included on the inside of the back cover of 
this report. 



The final part of this document treats specifically the thaw-consolidation 
problem which arises during the spring melting period and with which pavement 
weakening and possible failure are associated. The problem is still in the 
formulation stage. Field investigations of the thaw'-consolidation process and 
an hydraulic analog apparatus are needed at this time to advance further our 
knowledge of the problem. , ■ 


PREFACE 


, r This- report. represents the results of the first -stage of an investi- 
gation which has the following, ultimate aimss » : , .?■ r : 

? ,a# ,, Development of improved methods- of computing depth and rate 
of freezing and thawing in soils,, 

, b„ Determination of the relationship between the load-support- 


ing capacity of soils during thawing and the pertinent affecting variables 
(herein called the "thaw-consolidation” problem ) f and the development 


of improved engineering design criteria therefrom,, 

Active work on the program was initiated by the Arctic Construction 


and Frost Effects Laboratory in May 1952« At that time general outlines 
of the needed studies were prepared. These outlines anticipated that 
the studies would be made in two distinct phases s (a) Theoretical and 
analytical investigations*, to be carried out by contract*, (b) Practical 
application of the results*, through formulation of engineering design 
criteria, to be carried out by Arctic Construction and Frost Effects 
Laboratory personnel. The outline of the basic physical concept of the 
thaw-consolidation phenomenon*, as given in the outline of desired work*, was 
prepared by Mr* K, A, Line 11, Chief*, Arctic Construction and Frost Effects 
Laboratory, on basis of previous personal observation in the field,. 

After review of the initial Arctic Construction and Frost Effects 
Laboratory outlines*. Drs» Harl P„ Aldrich, Jr c and Henry M* Paynter of the 
Massachusetts Institute of Technology submitted a definite proposal 
presenting their proposed methods for making the desired theoretical and 
analytical studies, A contract was prepared*, and the studies reported 
herein were carried out during Fiscal Year 1953 » 


The report itself, which was prepared by Drs. Aldrich and Paynter, 
is self-explanatory of their work under the contract. During the course 
of the studies, the contractors held frequent conferences with Mr. James F, 
Haley, Assistant Chief, Arctic Construction and Frost Effects Laboratory, 
at which agreements were reached from time to time concerning the direction 
of the continuing studies. Under direction of Mr. Haley, the Arctic Con- 
struction and Frost Effects Laboratory also conducted laboratory tests under 
controlled conditions to check the validity of the frost penetration formulae 
developed by Drs. Aldrich and Paynter and made statistical studies compar- 
ing the theoretical frost depths with depths determined from field measure- 
ments during previous frost investigational programs. 

In the present report Drs. Aldrich and Paynter analyze in detail the 
penetration of freezing below the ground surface and present results of 
their initial studies of the thaw-consolidation problem. The investiga- 
tions will be continued by Drs. Aldrich end Paynter in Fiscal Year 195Us 
it is planned to make a more comprehensive attack upon the thaw- consolida- 
tion problem and to extend the present study of freezing to a consideration 
of the methods of computing depth and rate of thawing. 


• 4 
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y PART I V f. INTRQPUC TION V r £ -i 


1 - 01 . CONTRACT. - 

On 23 June 1953? the writers entered into a one year contract with 
the Corps of Engineers, New England Division for the frost studies presented 
in the following pages. This is the final report covering all phases of that 
contract. No. DA-19-0l6-ENG-23lUo 

The following paragraphs, which are quoted from Appendix "A” of the 
contract, describe the investigations which the writers agreed to conduct: 

1. ” The Contractor shall prepare nomographic charts for use 

in estimating the depth of frost penetration beneath paved and unpaved surfaces, 
and shall make a comprehensive theoretical analysis of the supporting capacity 
of thawing soils to determine the relationship between rate of thaw, consolida- 
tion and strength. 

2, During the first phase of the investigation, the Contractor 
shall review existing literature on mathematical methods of computing frost 
penetration with the objective of determining which of the various formulae is 
most suitable. The Contractor shall then prepare charts to permit the solution 
of this formula without the need of laborious computations and still include the 
fallowing variables affecting the depth of freezing, air freezing' index, mois- 
ture content of soil, thermal properties of soil (i.e., thermal conductivity, 
specific heat, and latent heat), freezing point of soil moisture and surface tem- 
perature transfer coefficient. 

3» The Contractor shall investigate more accurate methods of 
analyses of frost penetration problem with a view to numerical and machine solu- 
tions of the resulting equations, giving due considerations to the effect of such 




factors as surface insulation and soil stratification. During the seoond phase 
of the investigation the contractor shall, by means of theoretical mathematical 
analyses, attempt to develop practical equations and analyses to express the 
void ratio or saturated water content of any type of soil as a function of time 
during the frost melting period." 
l-oe. SCOPE. - 

While the proposed scope of this investigation is outlined under Sec- 
tion 1-01, in actuality the scope has been considerably broader. Principal 
additions covered in this report are (l) the presentation of a new formula for 
predicting the depth of frost penetration and (2) results of statistical stud- 
ies conducted to verify the formula. The SYNOPSIS serves to point out other 
factors considered within the scope of this investigation. 

1-03. PERSONNEL. - 

The contract studies and the preparation of this report were carried 
out primarily by the writers referred to as the Contractor in the contract. 
Several students at the Massachusetts Institute of Technology shared in con- 
ducting statistical analyses, performing computations, and preparing drawings 
for monthly progress reports and for the nomograph. 




•'!; - PART II, ANALYTICAL STUDIES ; :■ 

2-01. MATHEMATICAL FORMULATION. - 

a. Synopsis . - This article treats the mathematical basis for 
studies of thermal diffusion, depth of frost and thaw penetration,, as well 
as the seepage diffusion concepts required for analytical studies of the 
thaw-consolidation problem. 

After stating the first principles from which all developments 
begin, both exact differential equations and approximate difference equa- 
tions are given which embody these principles. Moreover, practical and 
fruitful fluid and electric analogies to the thermal problem are outlined. 

In order to overcome the handicap to analysis imposed by the com- 
plex latent heat conditions, several practical machine and hand methods of 
computation are presented. 


b. Fundamental Principles » - At any point and instant in a 


thermal medium, such as soil, the local rate of change of total heat con- 


tent (or internal energy) u must exactly balance the net efflux of heat. 
This amounts to a statement that heat is conserved and in the absence of 
heat sources, its flow must satisfy the continuity equation;? 


d u 

at 


div q 


0 


0 0 3 (2 — 1 ) 


Moreover, this directed heat flow q is experimentally found to 


be proportional to the negative temperature gradient, with the factor of * 


proportionality defined as the thermal conductivity k of the substance. 
This heat flow law thus yields the second fundamental equations 


3 


q 


«k grad v 


o o o o 


(2 - 2 ) 


where t is the instantaneous temperature field in the material* 

These two principles were first set forth by J« B 0 Fourier in 
his classic treatise on the "Theorie Analytique de Chaleur 1 * in 1822, and 
they form the basis of all rational investigations of heat conduction 
phenomena,, 


The heat content u of a substance is found experimentally to 
depend primarily upon its temperature v, and one finds under constant 
volume condi tions 5 at leasts that from a consideration of energy conserva- 
tion, it must vary only with the temperature* For example, the experimental 
curve of the internal energy u of water as a function of temperature v is 
shown as Figure 1 0 Three distinct, effectively linear segments are dis- 
tinguished, as marked by the numbers I, II and III 0 

The region I denotes the fluid state of water in which the tem- 
perature v increases approximately in proportion to the increase in stored 
heat u« This proportionality or slope we may define as the volumetric 
specific heat (under constant volume) C in the unfrozen or liquid phase, ors 

c. - 

u 3 v 

unfrozen 


a 


9 0 0 0 


(2 - 3 ) 


In the region 
phase and the 
specific heat 


marked ill, the substance (water) is in the frozen or solid 
slope of the curve may be measured by the frozen volumetric 
C^ 3 where s 



0 0 9 0 



h) 


C3 




- INTERNAL ENERGY DIAGRAM FOR WATER 

FIGURE I 


- 5 - 


> The flat or horizontal portion II of the internal energy curve dis 
plays the heat content -which must be released or absorbed to produce the 
phase change, and which is defined as the latent heat of fusion L of the ma- 
terial. 

As stated, the curves I and III are not exactly straight, but the 

deviations are usually negligible such that the small errors introduced by, 

assuming C and C both constant are consistent with the assumption of neg- 
u f /' 

ligible unit volume change, upon change of phase# 

Under these conditions, except in the region of the fusion (or 

freezing) temperature, the variation in u with time may be directly related 

to the variation in v, as the followings 

d u _ „ d v *#o#(2 - 5) 

d t & t 

where C becomes or C^, depending on whether v is above or below freezing. 


respectively. 


Substituting this last expression, together with Equation (2 - 2), 


into Equation (2 - l) yields 


C + div (»k grad v) = 0 

O t 


0 9 0 0 


(2 - 6 ) 




c# differential Equations # - If the substance is homogeneous 
such that C and k are constant -throughout a region, and is entirely either 
above or below the fusion temperature, then the following differential equa- 


tion holds for the temperature field v 

\ v 


dt 


- x 


V 2 ▼ 


«», »(2 - 7 ) 



This last expression may be re-arranged to give the thermal diffusion equa- 


tions 


v 

a t 


a V 2 


V 


» « » « (2 - 8 ) 


in which a 35 k/C is called the thermal diffusivity. For three-dimensional 


(x, y, z) problems p this equation -would be written 


d ▼ a / ^ 2 v 

at 




N 2 

d v ) 


d z 2 


« • • « ( 2 — 8a) 


For two-dimensional (x, y) situations it becomes 

a v _ 


a (~r~2 * 




O O 


00 (2 - 8b) 


while for one-dimensional (x) problems there remains 

a v _ a 2 v 


V _ 

ar " a 


a 


X 


o o o o 


(2 - 8c) 


do Approximate Equations «> - While the expressions in the pre- 
vious sections describe the frost penetration problem in a mathematically 
correct fashion, exact solutions can be found only for a small number of 

ts v 

idealized cases,, due to the complex conditions of the latent heat transfer 


and other effects <, 

When such conditions arise in engineering problems., it is the 
customary practice to resort to largely empirical techniques to obtain 
reliable prediction formulas<, However s a useful alternative follows from 
the substitution of an approximate system of simpler mathematical equations 
which may either be solved directly in closed form,, or which may be solved 
by various forms of computation or analogies 0 





In particular, it is often found effective to replace the con- 
tinuous Conduction equation (written in one-dimensional form for simplicity): 


q = 


a 


O X 


..,.(2 - 9 ) 


by a "lumped" approximation for a finite layer n of thickness^ x = ^ n j • 


in the form: 


- 


▼ \ 

( n+1 ” n) 




....( 2 - 10 ) 


or, by defining a thermal resistance R n = 


V k n 


ass 


9 . + (\ 

XL r — 


ill) 


....(2 - 11 ) 






The thermal continuity equation is here considered to hold for lumps of 
material concentrated at the terminals or end points of a series of such 
resistors; the relation between the net heat input and the stored heat for 


the n-th such lump would then be given by the equations 


du 


n 


dt 


■ V 


1 -*n 


o O o O 


(2 - 12 ) 


and the relation between the temperature v and heat u^, is indicated by the 


forms 


v n“V ( V 

where the function f n resembles the curve of Figure 1*. 


q o o 9 


(2 - 13) 


For hand computations and digital computing machines the continu- 
ous thermal continuity Equation (2 - 12) may be expressed discretely in 


terms of the finite difference approximation 


du 


u 


n 

dt 


n (t + At) 


© 9 0 • 


( 2 - Ik) 


A t 


to gives 


u n (t-Hlt)-u n (t) “At (t) - ^(tjj 

wQ®* 


4 0 9 


. (2 - 15 ) 


O 



In terms of a three-dimensional space lattice, the differential 
Equation (2- 8a) may be approximately represented by the finite difference 
equations 


At 


At 


a 


T l + V 2 


2v 


fax 2 


V + V - gv 

a 3 4 o 


V + V - Pv 

5-6 < 




(AZr J 


a © o o 


(2 - 16 ) 


■where s 

V At 

v(x 9 y,Z , t + At) 


It 

y,2, t 

) 


V 1 = 

v(x + A x, y,Z s 

t) 


v 2 * 

v(i-Ax, y»Z, 

t) 


V" 

v(x, y + Ay, H 5 

t) 


u 

v(x, y-Ay, Z , 

t) 


V = 

5 

v(x, y,H + AZ, 

t) 


1! 

v(x, y, H-A2 , 

t) 

If Ax 

= Ay 

= AH = A s, this 

_ ,a At.. 

' ▼ “ ( ) 

A 2 

A s 

may be abbreviated 

(T 1 + t 2 + T ? + V 


Y ° 

At 

or, with 

6 = a At/A s 2 , 
6 

to the expressions 


o o o o 


▼ - 6v ) 

6 o 

(2 - 17) 


V 

At 


pi 


v + ( 1-6 6 )v 
k ' o 


O O o o 


(2 - 17a) 


k=l 


In the same way. Equation (2“8b) becomes s 


1+ 


\t 


d 


\ + ( 1“U^)t 


k=l 


O 6 U 0 


(2 - 17 b) 



aad. the one dimensional Equation (2«8c) becomes i 


At 


= £2 


k=l 


T + (1-20 )v 


e • 


..(2 - 170) 


= (tj + v 2 ) + (l-2f)v 0 

t: has been demonstrated that certain stability limits impose conditions on 


the maximum values of ^ which can be used for computation (l) (2)** 

e, Fluid and Electric Analogies . - It is of significant interest 

\ 

to compare the thermal principles and equations to the analogous set of equa- 


ions governing the diffusion of water pressures in a confined fluid medium, 


This analogy gives rise to a simply realizeable physical model as described 

? . - ■ t 

1 

in Sections 2-01-f-(l) and l;-Oli-a below as well as in References 3 and J+« 
Moreover, the same equations hold for the diffusion of currents, and 
voltages in an electrical conducting medium which has a fixed capacitance 
to (ground. This, too, may be realized in simple one and two-dimensional 
passive electrical circuits, as described in References 5* 6 and 10 in ad- 
dition to the lumped equivalent as represented by the electronic analog 


computer outlined in Section 2-0L-f-(2). 

These analogies may best be summarized in tabular form as shoivn in 

Table I, The practical value of such analogies becomes obvious when the 

much greater manipulative flexibility of the fluid and electrical systems 

are considered. Moreover, instrumentation is significantly simpler and 

more precise in the analog media than in the prototype thermal system* 

However, one feature is distinctively absent, under normal conditions, in 

the fluid and electrical systems that indeed characterizes the present 

# 

thermal problem and makes studies difficult^ the presence of state changes 


♦Numbers in parentheses refer to Bibliography 
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TABLE I 


THERMAL - FLUID - ELECTRIC ANALOGIES 


ITEM 

MEDIUM 

THERMAL 

FLUID 

ELECTRIC 

A-Variables (l) 


Heat u 

Volume S 

Charge Q 

(2) 


Heat Flow q 

Flow Q 

1 

Heat, 0 " H 

Current i 

(3) 


Temperature v 

Voltage e 

B-Principles : 



fmm 


Continuity 

(i) 

+ v-q = o 


^-3 + V- i = 0 
<5 t 

Conductivity 

(2) 

q - -kV ▼ 

Q - -k^ H 

i = — (T Ve 

Capacitance 

(3) 

Aeu = Cdv 

dS = AdH 

dQ = Cde 


V 


11- 













and the concurrent latent heat effects „ The analogy to this situation 
■would require the equivalent capacitances to remain temporarily infinite 
until a corresponding amount of fluid volume or electrical charge had been 
transferred to or from the element* How this feature is secured in the 
hydraulic analog and in the electronic analog is described briefly in 

a , : .. , . . . . . , . ■ . , • ' • ‘ ' 

Sections 2-01-f«*(l) and 2“01«“f~(2), respectively* 

f • Practical Methods of Computation * - Due to the complexity of 
the latent heat conditions in addition to the normally heterogeneous nature 
of the soil, and the randomness in the surface temperature distributions, 
it is not possible to obtain very many practically significant solutions 
to these problems in closed, mathematical form* 

In order to overcome this handicap to rational analysis, several 
useful machine and hand computation methods are presented in the following 
section which are based on the principles outlined in the previous sections® 

(l) Hydraulic Model* - By applying the analogy relationships 
of Section 2-01-e and the ’’lumping’ 1 approximations of Section 2“01«-d, one 
finds that it is possible to solve a wide variety of thermal problems to an 
excellent approximation by means of a simple hydraulic model* 

Basically, the hydraulic analog would consist of a series of small 
vertical chambers or wells connected to each other at the bottom through 
laminar tubes or orifices* Each chamber represents the heat content, both 
volumetric and latent, of a lump or layer of soil at a particular depth* 

k ( r ■ - . ■ • . . . . , ... 

Distance along the model corresponds to depth in the soil* The surface 
level of the liquid in the well is analogous to the temperature* To. simu~ 
late the latent heat content of the layer an expansion chamber of suitable 



area would be provided at a level representing the freezing temperature of 
the soil moisture. The size of this expansion, as well as the size of the 

i 

vertical chambers, would be adjustable to correspond to the variations in 
latent heat and volumetric heat, respectively. More specifically, the 
sectional areas of the chambers are directly analogous to the C and L 
values, since storage volume S represents heat u and level H represents 
temperature v. 

The interconnecting capillaries or laminar orifices are 
directly analogous to the conductive paths with the fluid conductance rep- 
resenting the thermal conductance. These conductances oan be made adjust- 
able in a simple fashion so as to readily permit wide variations in the 
values of k between layers. 

The surface temperature variations are directly represented 
by variations in the level of an overflow-source tank at the end of the 
model which corresponds to the ground surface. 

The flexibility of such a scheme is directly evidentj; more- 
over, the same model can be used to represent consolidation, seepage flow 
and other diffusion problems. Such a model had been proposed several years 
ago by Barron (1+) and has been successfully used for latent heat problems 
by Backstrom (7). For diffusion problems associated with consolidation 
phenomena such a model has been effectively employed by one of the authors 
of this report (3)* 

(2) Electronic Analog Computor 0 « The hydraulic model just 
described becomes complex for 2-dimensional problems and 3 “dimensional 
representations seem entirely impractical, due principally to the difficulties 
of interconnecting the chambers in a manner which will still permit direct 

- 13 - 


observation. 


Jtt, 


These disadvantages are not present in an electronic com- 
putor operating on essentially analogous principles* A 2-dimensional 
array of elements can be mounted on relay racks covering a relatively small 
amount of wall space* 

The basic computing elements can be designed to solve the 
"lumped" approximate equations - the same as those which govern the be- 
havior of the hydraulic model* For the sake of simplicity,, only the 
circuits for a one-dimensional homogeneous computor will be outlined 
here, although, as stated, this method readily generalizes to the 2 and 
3-dimensional cases* 

‘ • fj 

In terms of operational blocks, the basic structure of 
the computor components representing one lump or layer of soil is shown 
in Figure 2A where the letters stand for the operations indicated beside 


the drawing. In particular, the "Z" - component represents the latent 


heat effect and has an input output characteristic similar to Figure 1* 

Thus each layer of soil is represented by an assemblage 
of computing components F which are interconnected as shown in Figure 2B 
to represent the entire soil mass* 

This particular type of electronic computor is referred 
to as an active representation in which all variables are represented by 
analogous voltages in order to expedite the instrumentation and inter- 
connection. Thus ground surface temperature fluctuations are represented 
by corresponding input voltage variations applied to the elements repre- 


senting the top layer of soil 0 
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Adding 

Component 

Output 
= Z(lnputs ) 

Coefficient 

Component 

Output 
= C * (Input) 

t 

Integrating 
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=/(lnput) dt 
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Figure 1 ) 

LAYER 
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v n + I 


II - II 


INTERCONNECTION OF F COMPONENTS 


SCHEMATIC BLOCK DIAGRAM 
ELECTRONIC ANALOG FROST COMPUTER 
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This particular representation is therefore a so-called 
"active” analog, in which all variables are voltages s in the absence of 
latent heat components, such computors are discussed in References 8 and 
9# This type of analog computor is distinct from the more common w passive B 
analogs (5) (6), in which temperature (v) is represented by a voltage (e), 
but heat flow (q), by a current (i)* 

Such a computor can be made to operate on so-called "fast- 
time” to permit display on a cathode ray oscilloscope, or it may be run 
at ”slow-time” (over a period of seconds or minutes) to permit recording 
on a single or multi-channel chart oscillograph* Actual surface tempera- 
ture observations may be traced by a stylus input device which drives the 
computor, permitting direct check against field temperature and frost - 
thaw depth measurements* 

The internal circuitry of a single "F” component is indi- 
cated in Figure 3® The input resistors (R ) and the feedback resistor 

i 

(R^) correspond to the thermal resistance* The feedback capacitor (C) 
corresponds to the volumetric heat of the soil layer* The dropping re- 
sistor (R^) determines the ratio of latent heat to volumetric heat* 

Since the computor may be arranged so that these circuit elements can be 
varied, it is possible to represent variations in soil thermal properties 
in a routine fashion* 

The only changes in the above circuit required for 2 and 
3-dimensional representation would be to provide four and six input resistors 
(R^) respectively* 


C»l6® 



R f 



CIRCUIT DIAGRAM 
SINGLE "F" COMPONENT 

ELECTRONIC ANALOG FROST COMPUTER 

FIGURE 3 



asj 


(3) I»BoM a Solutions* - As a means of obtaining reliable 


rational data on the freezing and thawing process as part of the second 


phase of the contract, and due to the non-availability of the hydraulic 
and electronic analog equipment just referred to, a number of computed 
studies were undertaken using the ’I*B*M* Card Programmed Calculator of 
the 1LI..T* Statistical Services Division* 

These studies had three main purposes, as follows, 8 

(1) Verification of the difference equations as outlined 
in Section 2-01-d, including the effects of latent heat* 

(2) Determination of the effect of the shape of the sur- 
face temperature-time curves on the total depth of freeze* 

(3) Investigation of the shapes of the depth- time curves 
during both the freezing and thawing period* 

The basic numerical solution techniques have been exten- 
sively outlined in the available literature (ll) (12) ( 13 ) (ll|.) following 
in all essentials the procedures outlined in Section 2-01-d* However, 
to the best of the authors v knowledge no numerical solutions have been 
published involving the diffusion process when latent heat effects are 
significant* 


In order to ascertain the precision of the numerical ap- 
proximation procedure it seemed desirable to run check solutions to com- 
pare with the rational analysis of Section 2-02-c* Accordingly, twelve 


\.r .. ... . , * . . ; *•,. •. ... - , 

test solutions with various values of 0( and yUu were carried out on the 


T*B.M» equipment* These are presented as Tables B-I through B-XII, in- 

# * t 

elusive, in Appendix 1 B, and assume a step change in surface temperature* 



Two solutions given as Tables B-XIII and B-XIV in Appendix B s dis 


play the effects of a sinusoidally trying surface temperature 0 They 
refer to purposes listed as items (2) and (3) above and are considered in 

i 

Section 2-02-io 

These fourteen solutions,, as •well as several other check 
and trial solutions not presen ted 9 were obtained by automatic computation 


of the following set of equations and conditions^ 


k+1 


* S 


n , k 


0 ( 


where 


e 


n, k 


9 n e k 


n. k 


* S 


n, k 


e , + e 

n-l, k n+1. k 


for S , > 1 
n„ k ^ 


2s 


a, k 


0 o 


us (2 *-18) 


* 0 


= S 


n, k 


for 0 


for S 


n s k 


0 


1 




* J> o • 


(2 - 19) 


The set of conditions (2-19) defines a normalised version 
of the internal energy - temperature relation similar to Figure 1„ with e 
being a measure of the temperature v and S a measure of the heat content u 0 
Rendering the relation in this simple form permitted a satisfactory pro- 
gramming of the calculation* However,, it makes necessary the proper scal- 
ing of variables and interpretation of re suits 5 as indicated below and in 


Appendix B, 


The dimensionless parameter ^ is defined in Section 2-01-d 


as 


0 


A. t 


A s 


where a is the ( thermal )diffusivity s A the time increment., and A. s 


S 


the depth increment. For stable calculations in the one-dimensional case ^ 
must be taken equal to or less than 0,5* Ibe values chosen for the various 
I,B,M* solutions are indicated in Appendix B, 

“ 19 - 



The relation Between the computed e - values and the actual 


temperatures v (t) is given by the expressions 


v (t) = v • • 

; ' ■ s v. 


e. 


k 




F e 


....(2 - 20 ) 


where v ■ F/T is the actual or equivalent step change in surface tempera- 

’*■ ■■ l, S 7 • 

; » 

ture below freezing and = Cv /L« In other words, the e - values 

. .... . , ■■ «... . ...■■ - . S 

bear the same relation to the parameter that the v - values bear to 

• ••»•• . * * . ■ 

v * The time t corresponding to any k - value is given by the defining 

5 


relation : 


t * kAt 


9 9 e 


(2 - 21 ) 


while the depth x corresponding to any n - value is given by the definition 


x = n/i s 


9 9 0 0 


(2 - 22) 


Thus instantaneous temperature gradients as well as tempera- 
ture-time curves can be plotted from the values given in Tables B-I to 
V-XIV just as for actual thermocouple data 0 Such curves were drawn by the 
authors and the conclusions are given in Section 2-02-f « 

(10 Hand Computations * > Using procedures very similar to 
those outlined in the previous section, it is possible to carry out numeri- 
cal solutions using a slide rule or desk calculator* In fact, for certain 
cases completely graphical procedures are practical* 

Efficient tabular forms for such numerical solutions are 
detailed in the available literature (ll) (12) (13) (llj.) while the Binder 
graphical procedure for heat flow problems is outlined by Jakob (15 ) 0 

However, the procedures of the cited references must be 


mod' 

the 


Lfied to account for the latent heat effects. This involves use of 
e - S relationship outlined in the previous section, A portion of a 


sample calculation sheet is given in Table II, The calculations follow 
dirjectly from Equations (2 — 18) and (2 - 19) as indicated* 




TABLE II 

SAMPLE HAND COMPUTATION OF FROST PENETRATION 
(A) NUMERICAL CALCULATIONS? 



0 

1 

2 

3 


0 

-loUo 

X 

-1.1)6 

X 

-1.50 

X 

-1.1*8 

X 


-l.o6 

- J 4 2 

- .90 

- .36 

- .79 

“ .32 

- «73 

- .29 


- .3 h 

- .314. 

- .58 

- .58 

- .71 

- .71 

- .75 

- .75 

1 

- .i;6 

- .18 

- .58 

- .23 

- .71 

*** 0 28 

- »75 

" .30 

, 

+ .12 

+1.12 

0 

+ o97 

0 

+ .81 

0 

+ .58 


- .08 

- .03 

- .18 

- .07 

“* .12 

— .05 

- 0O9 

- .01* 

1 

+ .20 

+1.20 

+ .18 

+1.18 

+ .12 

+1.12 

+ .09 

+1.09 

3 

- .02 

- .01 

- .02 

- .01 

- .06 

- e 02 

- .05 

- .02 


+ .22 


+ .20 

X 

+ .18 

X 

+ aU 

X 


(B) KEY TO CALCULATIONS s 


n-1 


n 


n+1 


k k+1 
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2-02. DEPTH OF FROST .PENETRATION# 




a. , r Syno psi3. - This arti cle outlines,, the wri t er si development of 


' V J' V' ■* ' '-? 1 




g;-:. JIi-> tW -Vn* . 1* ’ !• w >i<- / >■ -Z I * 


a rational formula to predict the depth of frost penetration which takes into 




c.\. rD. { 

c. t-i* • * “+> v* . t\ 


consideration nearly all of the factors of primary significance to the prob- 
lem. The nomograph included in this.. investigation employs- this formula as 


its" basis . 


f r* v* J 


7** 

k i . . 

.< » ** 

V » •J! /. «*.. • ’ 


■, This formula is premised on obtaining the exact solution to an 

V ■ j* ‘ 

'"if- • v- s*-- •• >»vf —• — _ * 

J* «*■ ‘ vrZi . 

idealized problem which is sufficiently representative of actual physical 

* \ 

conditions to serve as a basic reference for any future correction terms. 

^ , Following the derivation of .. the xationail;fbrmula a the writers * 
present the results of their statistical studies of actual frost penetration 
data. , It ..is conclude4,>that the .rational .f ormula generally yields better re- 
sults .than, otheri, expressions.., which were, .studied... .... 

.... The. latter, Parts, of y this .section, are devoted to- a ,. study of,, the, 

« 

effect of a cyclic surface temperature -fluctuation on the . depth of frost ..... 
penetration., .and, the, rate.,,, of freezing.. and,thawing«>,, r .Finally, - .the multilayered 


system is treated 


... % ar* .cs .Si>: jxc 5 i:o cr.ii : .-■o.o.ar \.c- .: x-< 


fr»- 7T * O i-i f , t 


b.* Review of r . Literature , --,As.,a. first - f s,tep;. in. these,. studies. 


existing literature pertaining to. T prediction- of , the 1 depth of frost penetra-, 

/ \ “* ' 

tion was reviev/ed*. .. p The. principle, .objective was. ..to determine the most suit- 
able, formula on which to base the nomograph, -Consideration.: was .given .to three 


criteria: 


' r ‘ 1 . ;■ r'’ i •• . v * . r 

’ r.‘ ? ; *2 ■’< > {-} •". r •** -t. i t ... . ,r ? •• i- ,t % 

^ O .v. v ‘J.* i O •/ 4 ,v “ .iA-V < 


, (l); The. mathematical validity of, the, existing formulas ; 
,, . i C2k--v?^ 1 ©a. e l ls l9i of ..application of. each., formulaic . 


iit s.dioeb six Ii..f!dm3v»*8.d;;e .$> iibMw 




I 


(3) The statistical reliability in predicting depths of 
frost penetration as compared to field observations. 

In particular, four formulas were considered, namely s 


X 


l±8kF 


0 0 


..(2 - 23 ) 


X = 


14.8k F 


L + 8d\&(C g + OoOlw) 


X = 


\! 


i*8kF 


L + C (v + — ) 
o 2t 


o 0 


• »(2 - 214 ) 


O 0 


• * (2 - 25) 


X = 


P)|kF 


L + C /v + F 
v o 


2F 


0 O O © 


(2 - 26) 


The first equation is the Stefan equation which considers latent 
heat only. The second equation is the so-called "Minneapolis” formula. The 
third and fourth expressions are formulas which have been studied by the 
Fhost Effects Laboratory of the New England Division of the Corps of Engin- 
eers, These last three consider certain of the significant volumetric heat 
effects by making various rough approximations. 

After a preliminary analytical and statistical review of the above 
four formulas, the writers concluded that Equation (2 - 25) was the most 


sound. Other statistical studies 


( 16 ) 


comparing actual depths of frost pene- 


tration against those computed from the above equations had indicated that 
Equation (2 - 25) gave reasonably reliable results. 

However, the writers were not entirely satisfied with this equa- 
tion, since it did not account for the volumetric heat effect in the unfrozen 
ground, which causes a substantial reduction in frost depths in the more 
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temperate jcliinat.es *, .Accordingly, , they undertook the, derivation <of a more r <1? 

rational formula as outlined in Section 2-02-c „and. Appendix A* ..^Thist was ,, 

do cj» ,*.C'L;rs£&:oc ro.oc -? aert.sw.i 

f 17) 

based on the. earlier .work of Neumann, as, reported in ..Cars.law, and . Jaeger V , 

dealing with melting and. solidification problems „in ,still mater . 

JLA*Oft i°s^ w Uk' ui x'i.l'i c\fter. .;o ©;>&'* xb -? xo eno ca Juiiip* -•-■• 


sVi* , . -•* „*• <• !' ■«* 

» . ? -5 ; .- f-V. '•, ■< W X> 


After this formula was developed, it was discovered to be. identi- 

ju.ox edcr xvxx; xx :x >:>.<¥ x.ss&r_ _ci_ja ei5t< •- vX 

(18) 

cal in nature to a formula developed by W* P„ Berggren' , which has, been ,,, 

n 

expressed in the forms. , > ... e , ...,>•*> 

- SB \J at 


f 


X ® 2B \/ at — ~ v ....(2 - 27) 

•j toia;.'. ; ; . ,</ 7 

by Shannon 7 • The relation between the factor B v and the notation of the 

t- *» 

writers is outlined at the end of .Appendix A* The unfavorable report by 

si* »’.i ‘xo >.^c^b T:-t .u<i-xys.^o; s jot 17 

Shannon in the reference cited had caused the writers to be initially, mi s- 

oo.'ix-'-v/-: 1 'x-.jO'd hJo"/.. ’ivv t xxx'aSv o.cxr _ o svi- '- O's/r;" 

- 'f ' ’ k 

lead; the 50^ discrepancies mentioned, in. connection with Shannon ? s calcula«. 

xsii'ri- -ixf; xxx x-' ; :r '.;:xxx?o-; '«;. li , »«o/.'ser? 

' ' " ¥ '~ *4,$ N *V* **■ f, ‘ W 

tions must be attributed largely to inadequate data, pn the. thermal, properties 


of the soils considered 




The statistical studies, undertaken by the writers and by the, 5 , 

~ -• V-' - - j ' 1 - 1 »*■..■. 

Frost Effects Laboratory as reported in Section 2-02-d~ £2) «;* testify strongly 

*•-•' S;-. to. ■ ,X 

to the general reliability pf this modified. Berggren formula* f . r 

0&^kP l ‘Ly’+ (’ r^S V; ' & \ ‘ *;vf V ( *CW .L- .. ; 

c. . Derivation and Interpretation of Rational' Formula, , - A de^ 


K' . X y.4 '•»* ^ • K. 


tailed derivation for the writers’ rational formula, for the depth of frost 

•xr.xox s ., -re- ! ?x fx.tutvt. exi.r yfnem si: d~Ts&:< -jin*# y cx : 

penetration is given in Appendix A «* The most useful form for using the, 

t’v&OU ^ '.O.'XS; etc?'. : Jlo L* <K ; ~£cvS 'IV T •* 

formula is given below as Equation ;( 2 - 29)* , , . 

a.i 0*1 <$r .•.v:,f«. xxrxhxx x;: 'X':: v';-: xxx x 1 x.xxx ; .ui xv 2 

As indicated in Appendix A, it is assumed that the soil is 

fxtx.xx:': ,» Zd Vv. ; rd y/ X J" : -~ 'fX:- 


V 


a semi- 

»• #*j. ■:> j'f'r?" 

« 4 4*1 O. i*- • . J 4 ‘J . 


infinite mass of uniform properties and having initially a uniform temperature* 

r Q-'X : y AlI’B /,' ■’> J .SC 

It is further assumed that the surface temperature, is suddenly changed from 

£ Qt •- HI* •!-- C f r. • " • <rw "' !r ^ ' te r '5.* liS«- 0-12' 

* . s> 

its initial value v above freezing to a temperature v below freezing* 

o s 


fH * 
•'* i- 


,K’a 
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The equations which must hold for this problem are the diffusion 
equations (2 - 8 c) in both the frozen and unfrozen soil as indicated in 
Figure I4* The latent heat property becomes a continuity condition which 
states that the rate of flow of heat in the frozen soil at the frost inter- 
face must be equal to the sum of the rate of heat flow in the unfrozen soil 
at the interface and the heat given off at the interface when the soil 
moisture freezes* 

The solution of these equations may be written in the forms 
\ , / UfScv t 

x ° Ay- 8 =.,.(2-28) 

For practical computations for the depth of frost it is generally 

assumed that v represents the average surface temperature below freezing 

s 

during the freezing period, Thus, if F represents the air freezing index 

and nF the surface freezing index, thens 

nF * v t 
s 

and therefore s 

...,( 2 - 29 ) 

which mil be called the modified Berggren formula* 

The term X , a correction coefficient modifying the square- 
root term, which is merely the familiar Stefan equation, is a complicated 
function of three dimensionless parameters 0 ( , yCL, and § * It has been 

shown that for all practical purposes § may be assumed equal to 1*0 in 
which case X is given graphically by the curves in Figure 5 as a function 
of 0 ( and where 8 

v t v 

thermal ratio, 0 ( - — 2 — = ° *, 0 *(2 - 30) 

nF nv g 

CqF 

and fusion parameter, /U* = — r— **o *(2 - 31) 
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The thermal ratio o( is a measure of the* ratio, of, the- initial.-,* 
ground temperature, or in practice, the mean annual temperature above the 

** • ! : A' *?- 

freezing point, to the time-average surface temperature below freezing during 

j. ' i 

. — ? f*'* ^ ■ /> r ‘ 

theCfreezingv period* It is apparent, then,- that - decreases in the colder 

i <• • /• i '*■ •%. f " $ 

\ ?“ : .. v / 

climates* In parts of Alaska Of may be zero or actually negative* 

■ V) V . r> 

; .i;r , The fusion parameter ZU- is* a measure '-bf- -the- heat removed in the 

’’ 1 f V * : V -V* :> 1/ 

^ V i*, .* * ^ t % *S 

frozen soil below the freezing point as compared to the latent heat of the 


soil moisture* As the latent heat tern becomes very large, /U* approaches 


.*ti ft .',* ' : , 


zero* It is seen, furthermore, that yCc, varies with the freezing index 


and is therefore largest in the arctic climates 
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d. Comparison of Rational Formula with Other Formulas* - A com- " 


pari son of Equation (2 - 29) with other formulas proposed and' used' to com- 
pute the dephh of frost penetration will be made ah two levels* First,' the 

.» / ' ’•? t'rf ,♦••• '*• />■ '**'»'* ' ;<;/ /C’i ‘fjy • 4 t;.. : -:*j /-v? ,5 7 'h. *' ’j 

assumptions made in each derivation are outlined and a general evaluation of 
the effect of these assumptions on computed results is presented* ' Second, 


results of statistical studies using actual frost penetration data are J 


analyzed. 




(1) General. - Equations (2 - 23) (2 - i25), and ^2 -' J 2§) a of'' 


i ’ , jJ ‘ ’ v ‘ '.I' r“. * *" s * v - f C ,1 _ 1 ' / 

Section 2-02-b will be used for this comparison. It can be shown that all 


■ft “ . 


. - » * - r . . .... A . t , -• ‘ ' y : **'■ *■ y ’ .i ' iT** * 

three expressions can be rearranged and written in terms of a correction co- 


efficient which is some function of o( and 
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previous seotiorLo Thast 


'X: 


«‘t, J ■/ ■** *' 5 " i y J ‘C x r 'i'- V ft .... 

as defined in the 
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3 -i-v / - i % 1 Ij_. /„ • / 1 . 
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if 



where X is given by the following expression for formulas (2 ~ 23) 3 (2-25) 


and (2 - 26) respectively: 

Ai “ 1,0 
A, = 


9 0 0 0 


(2 - 32) 


VTT + 0.5) 


OOOO 


(2 - 33) 


A 


0.707 


\[ 


0 0 0 9 


(2 - 3U) 


1 + yU, ( C< + 0.5) 

Curves for ^Vi* X 2 and X^ are shown in Figure 6B* 

Since the Stefan equation (2 - 23) ignores the volumetric heat 
given off as the soil temperature is lowered to and below the freezing 
point, is equal to zero and therefore. A. = loO for all C>( * 

Computed depths of frost penetration based on this formula will generally be 
too great a This is especially true in temperate regions where actual values 
.of A may be 0*6 or smaller* Prom Figure 6 a it is seen that in temperate 
Kansas, for example, the correction coefficient is perhaps 0*55 ‘Which means 
that the actual depth of frost penetration would be about 55$ of the value 
determined from the Stefan equation* On the other hand, in Alaska the equa- 
tion would be expected to yield good results* 

0 

Equations (2 - 25) and (2 - 26) consider the effects of volumetric 

- . * z' 

heat under various assumptions* They should be expected to yield somewhat 
better predictions than the Stefan equation* It can be seen by comparing 
Figures 6A and 6B that for small 0( equation (2 - 25) gives good results 
while equation (2 - 26) with a much smaller A will predict too shallow^ 

An engineer in Kansas or even Nebraska where 0( is large would find that 
equation (2 - 26) was very good while equation (2-25) predicted too deep* 
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The statistical,, studies described in th e H f ol lowing .section 0 generally 


confirm the, above statements* 
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(2) Statistical Studies* - r The following studies ~ are aimed,, 
at comparin g pr edicts d depths. .of frost .penetration gi yen - by .the .rnodif i ed 
Berggren .formula (2 - 29) with actual -depths. . determined ;; bjr .themocpuple read- 
ings and test pits* ,. These results are compared .with, predicted re suits, r .nsing 

‘ ’ : \ T — ‘ , * “ " - -S 

Equations, (2 -25) and. (2 “, 26 ), the latter including the rpaye^e$b ^l^clnqe^Sifl* 
Data for this. study- have been taken from 1±8 cases given in reference, (16.)* 

Two studies were conducted by the writers* Average values,., of C* i ( L 

; XHs 17./ y * -J v ^ >*i .'.J .• « ;* ■ , .... ; . 4 V : r-c •* ^ rov..,^ ♦'•uv 

and^ k_as given in ref e real ce _ ( 16)_ .have bean . used* ^_addit.lon,_..^ 

ing .index .has been, used for both.. s.tu dies* .. — ] 

£ ? '.“'T-.-' ■V.'vV^ A ^ | 

The first analysis, involving the modified Berggren formula only, = 
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The graphical representation of the dispersion is shown in Figure 7» 
This may be compared with plate 5 of reference (l6)* 

In general these results do not show a prediction which is signifi- 
cantly superior to the other equations® 

The. second analysis is based on inches variation between observed 
and predicted depths 9 comparing Equation (2 - 25 ) 3 Equation (2 - 26) (with 
pavement thickness added) and Equation (2 « 29 )o For this study the data 
have been arranged in two groups 9 one with the 2 1\. highest values of the 
actual correction coefficient a * (<>915 to *715) and the other with 

the lowest 2 1+ (o7li+ to ®295)<> The following results w ere obtained s 



Group I 

/ 

Group II 


X s (.915 to .715) 

(cT'lij. to « 

295) 


Avg® Obs<> Depth = 

52®3 inches 

Avgo Obs© Depth ~ 

l;0o9 inches 


Avg® Deviation 

Standard Dev® 

Acg® Dev® from 

Std c Dev® 


from Observed Depth 

from the AvG® 

Obsrvd® Depth 

frm c the Avg® 

Eq, (2 - 25) 

+ 5 . 4 " 


+19o2 n 

^9.3" 

Eq. (2 - 26) 

- 4 . 4 " 

!5o9 w 

+ 7.1" 

+8 

Eq. (2 - 29) 

- 0 . 4 " 

+1*5" 

+12® 2 n 

+8ol n 


It is apparent from the above table that- in general for high 
values Equation (2 - 25) predicts frost .depths which are too great while 
Equation (2 ~ 26) predicts too shallow*, For low t all equations pre- 

dict too deep® Equation (2-25) predicts about 50$ too deep while Equation 
(2 - 26) is less than 20$ in error® These .results* then,, confirm the general 
comparisons made in the previous section® 


A 


observed depth of freezing 


y IjBkF 
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PERCENTAGE DEVIATION FROM OBSERVED VALUES 


DISPERSION OF COMPUTED FROST DEPTHS 
USING MODIFIED BERGGREN FORMULA 

FIGURE 7 
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t Tq be more specific, one may readily see from the above table, r r .; f - 
that for the first 24 points ( Group I ) the modified Berggren formula (2 - 28) 
predicts the depth of penetration within an average deviation of ; -0.1+ inches 
out of a group average observed depth of 52 inches, with a probable error 
of +3.0 inches (0*6714-5 x 4*5)* Thus the average deviation for the data of 
Group I is less than 1 % of the average depth and is better than 10 times 
: smaller for the new formula than for Equations (2 - 25) or (2 - 26). This 
means that, on the average, for values of X between 1.0 and 0 ,707 ((say), 
one may expect Equation (2 - 29) to predict about 82 per cent of the obser- 
vations within 6 inches of the actual values. For example, in Group I of 
the data studied, 21 out of 24 points or 87 per cent were within 6 inches.. 
For the second 24 points (Group II), Equation (2-29) predicts 
the depth of penetration within an average deviation of +12.2 inches out of 
a group average observed depth of 4l inches with a probable error of *5«5 • ; 
inches. Then, on the average, for values of less than 0<>707 a one may 
expect the new formula to predict only about 22 per cent of the observations 
within 6 inches of the observed values. For example, in Group II, 5 out of 
24 points or 21 per cent were within 6 inches. v.v . -,v -x- ed 

....... Since, in Group II, for all the formulas, both the average devia- 
tions and the probable errors are significantly large, one is led to sus- - . 
pect that at least part of the fault lies with the assumed data. From a - v , 
superficial analysis of this group, it would appear that the actual values 
of the thermal conductivity k are in some cases significantly different 
from the value assumed for this and the previous studies reported in 
Reference (16). Moreover, the assumption of k u 88 k^, seems to be appreciably 
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in error for certain of these points 8 These conclusions would suggest 
that for Those cases in which the approximate values of CXC and ytX* point 
to values of less than 0»707* & more refined determination/ of reliable 

soil constants is in order 0 


Following a suggestion by the writers ,, the Arctic Construction and 
Frost Effects Laboratory made new statistical studies using for each case an 
average thermal conductivity k determined from the water content and dry 
density of the soilo Results of these studies showed that Equation (2 - 29) 
gave predicted results superior to other formulas 0 The Arctic Construction 
and Frost Effects Laboratory therefore approved the writers* r e commend at i on 
that the nomograph called for in the first phase of this research be based 
on the rational formula (2 ~ 29 ) developed by the writers s 

So Limitations and Applicability of FormulaS o « Before discussing 
the limitations of Equation (2 - 29) it is of value tc list again the basic 
assumptions made in the derivation of the formula z 


(1) The soil mass was considered to be homogeneous and one 


dimensional 


(2) The initial temperature in the soil mass wa« assumed to 


be everywhere constant at the value v above freezing,, 

c 

(3) The surface temperature was assumed to change suddenly 

from v_ above freezing to v below freezing,, and remain steady at this latter 
o s 

value o 


(1*) The effect of latent heat was considered to introduce a 
heat source (or sink) at the moving freezing point interface., with a heat 
flow proportional to the speed of motion 0 
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(5) Under the above conditions an exact solution was ob- 
tained, taking into account both volumetric heat and latent heat, as well as 
the variations in thermal properties upon freezing of the soil moisture* : 
This exact formulation of the above idealized case was expressed 
in terms of three parameters,, namely the thermal ratio,£fV the fusion 
parameter,^^. , and the root diffusivity ratio, & « These same con- 

stants would govern the behavior in soil freezing problems under the much 
more complex conditions encountered in actual practice* Thus it may be . / 
expected that the depth of frost penetration in these actual cases is af-^ 
fected in much the same way by Of , , and <§" , as the idealized 

case for which the formula was derived. 

However, the failure to realize these simplifications may produce 
significant divergences; therefore, the most important of these deviations 
are outlined in the following paragraphs* 

The actual surface temperature variations are not uniform over the 
surface area, but rather are appreciably affected by ground cover, pavement 

type, etc*, as well as by local variations in wind velocity, shielding, 

' ■) 

snow cover, etc. For any measure of success using one-dimensional methods. 

of prediction, considerable care must be taken to derive a set of consistent 

and truly representative conditions and equal care must be taken in the inter 

pretation of formula results for actual cases, 

All actual soils are layered or otherwise variegated to a greater 

or lesser extent and the assumption of homogeneity must be accepted a^-au^ 

mere approximation. As the variations from layer to layer become; 

pronounced the formula predictions may become seriously in error unless 
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deliberate care is taken to compensate for the layering as outlined in Sec® 
tions 2-02-g and 3-02-c , 

The initial temperature in the soil mass upon the start of frees* 

ing is not uniform at a value v above freezing, but rather varies in some 

o 

non-uniform way from the surface temperature down to a temperature comparable 

k 

to the mean annual temperature at a point sufficiently deep in the soil. Rep® 

resentative of such variations are the values at the start of freezing given 

in Tables B-XIII and B-XIV, However, as testified by the I. B*M. solutions 

themselves, as well as field data and the analyses of Section 2-02-f, the 

effects of such non-uniformity of temperature on the correction factor, , 

are not appreciable so long as v is taken as the mean annual temperature. 

o 

Much of the same reasoning applies for the shape of the surface 

temperature-time curve and its effects on penetration depths, as discussed 

in Section 2-02-f. Again, as long as the equivalent step change v is taken 

s 

as the mean temperature below freezing during the freezing period, a good 
agreement exists between actual seasonal depths and those predicted by the 
formula. 

Other limitations and modifications of the formula as applied to 
practical problems of predictien are discussed in connection with the nomo* 
graph in Sections 3-02 and 3 - 03 • 

f. Depth -Time Curves . - "While Equation (2-29) gives consistent 
predictions of seasonal depths of frost penetration, attempts to calculate 
the depth-time curves based on partial freezing indices are likely to en* 
counter substantial errors. This may be demonstrated roughly as in the 
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following paragraphs. 



If the soil is initially at a uniform temperature y and a sur- 


face temperature variation v (t) below freezing is applied, then to a crude 

s • ’• ■* . ‘ • . 


V ** \ ‘.t 


first approximation the following equation may be shown to hold: 



....(2 - 35) 


This equation accounts for the latent heat effect, and the most 

It may be normalized to 


significant terms due to volumetric heat effects, 
the form: 


(t) - -iM. . 

Jr ' 


h l +W 


m 



m 

2 



1*8 KF 
L 


r< 


i+pe 


■)d7r 


i+cr 


i +pe 


) its' 


...=(2 - 36) 


• • 


..(2 - 37) 




(2 - 38) 


V 


where CP = = v»A 

P /U /!. ■+*/*■- 


6fr) 


v /(P/T) 
s 


T = Duration of Freezing Index (in DAYS) 

T 

F - f v (t) dt * Freezing Index (in DEGREE-DAYS) 


J 


s 


* . ; ^ ■ t/T « Relative Time during Freezing Period 

(i*e* 7^® 8 corresponds to start of freeze 

"■ I corresponds to end of freeze.) 

Then jzf ft) * (Zri/j m measures the relative SHAPE of the depth-time curve 


From the definitions above it is necessary that: 

'ec/2r = O 

'0 

w^8 a 


I 



whence it becomes evident that the function £ C s also, merely 


measures the SHAPE of the surface temperature curve during the freezing 
period* Thus we may investigate simply using Equations (2 « 3&) and 
(2 - 38) the effects on the depths time curves of various surface tempera- 
ture time curves* 

The, simplest such case is that of the constant step changes* 

In this case £ * 1,0 * constant throughout the period and the in- 


tegral for 


becomes 


V 2_ I f' jTr 

| + 6- J a I + 


Cl + o%i + P) 


0 0 9 


*(2 - 39 ) 


Upon substituting the defining constants for the parameters O’ 
and p s this expression for A is seen to give the value corresponding 
to Equation (2 - 25 ), namely s 


X- 


\] I +^C^+o.5) 


o » .o -o 


(2 - I4O) 


Thus the problem of investigating the effects of the surface 
temperature curve on the frost penetration depth curve may be factored into 


two parts, namelys 

(a) The SHAPE effect 

(b) The MAGNITUDE effect (A) 


Equation ****(2 — 3&) 
Equation o**o(2 - 38) 


It is easy to show that for any shape €&) symmetric about 

the midwinter point 7)* = 0*5$ which is very nearly true for all practical 
instances, then may be very closely related to the mean value £ • l o 0 


by the following expressions 

A 2 ^ 1 


i 


I +<3- 


i + pe 

■ 39 - 


O+^Xi+fb 


O O 0 6 


(2 - i| 1) 



This fact would indicate, that the shape of ; the surface temperature 
curve during the freezing season has, only a very slight effect on X , and 
therefore on the total depth cf freeze© This hypothesis has been substantiated 
both from field measurements and by the studies referred to in Section 

2-01-f-(3)© 

On the other hand,, a slight effect on the shape of the depth-time 
curve does. exist and may be explored by substituting the above expression 
for ^ into the equation for (£> to give the forms 


o o e <> 


(2 - 1 * 2 ) 


The integral may be approximated by a finite sum in the forms 


0 *- 


O O 9 * 


(2 - 1 * 3 ) 


where h = and where jf'ediy — |.0 is replaced by the 

condition | q * Practical values of p are of 

magnitudes not much greater than 0.2© Accordingly, the table below has been 


calculated with h =s O. I and values of £ and as indicated, in 

order to show the effect of P on (£) © 

. From Table III, it is clear that the parameter p has no sig- 
nificant effect upon the shape of the depth-time curve over its physically 
practical range. This permits one to draw valid conclusions as to the shape 
of the depth-time curves from the simple case where only latent heat is im- 
portant Since 0^ is then merely proportional to the freezing 

index, such cases may be explored very quickly© Conclusions based upon this 

approximation are shown in Figure 8© 
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TABLE III 


SHAPE OF DEPTH- TIME CURVES 
VALUES OF 

BASED UPON A PARABOLIC DISTRIBUTION OF 


k 


VALUES OF P 

0 

0.1 

0.2 

0 

0 

fr 

1 

0.29 

0.17 

0.17 

0.18 

0.19 

2 

O.76 

0.32 

0.32 

0.33 

0.35 

3 

1.12 

o.ltf 

0.1*7 

0.1*8 

0.1*8 

u 

1.36 

0.59 

0.60 

o.6o 

0.60 

5 

14 a 

0.71 

0.71 

0.71 

0.70 

6 

1.147 

0.80 

0.81 

0.80 

0.80 

7 

1.36 

0.88 

0.89 

0 o 88 

0.88 

8 

1.12 

0.95 

0.9l* 

0.9l* 

0.91* 

9 

0.76 

0.99 

0.99 

0.98 

0,98 

10 

0.29 

1.00 

1.00 

1.00 

1.00 
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LINEAR 

► TIME 
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STEP CHANGE 


TIME 


THAWING CURVE 
APPROXIMATELY 
LINEAR 


FREEZING CURVES : 


(a) FOR PARABOLIC TEMPERATURE: 

APPROXIMATELY PARABOLIC WITH 
VERTICAL AXIS THROUGH B 3 


(b) FOR STEP TEMPERATURE: 

APPROXIMATELY PARABOLIC WITH 
HORIZONTAL AXIS THROUGH A g 


SHAPE OF DEPTH-TIME CURVES 

FIGURE 8 
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.. g.. Multilayered Systems, - In any. given layer, m, the instantane- 
ous heat flow, assuming a linear gradient, is given by Equation (2 - 10) or: 

A v m .. ..... : ... .. . • • 

•• m ,V.; :• • • . • .. 


q = -k 
4 m m 


*m 


• • • • (2 “ 1 + 1 +) 


which may be rearranged to solve for the temperature difference in terms of 
flow and the thermal resistance of the layer, = d , to give: 


A ▼ * R q 

m mm 


• •••(2 - 1 + 5 ) 


For an entire series of layers, from the surface to layer n, the total tem- 


perature differential may be given by: 

n 


n 


( 




R ) 

nr 


• ••(2 - 1 + 6 ) 


If the n - th layer is at the freezing temperature, then v n = 0, 

Moreover, if the only source of heat is considered to be the cooling of the 

layer from v to the freezing point and subsequent revival of the latent 

o 

heat L n , then the total heat removed to freeze the layer n is given by: 

U n “ j\ + Vo] d « ....( 2 - 47 ) 

Thus the equation for the freezing of this layer can be written: 
du 

q a — XL- « • . • ( 2 “ 1 + 6 ) 

^ dt 


or 


dU = q^ dt 
n n n 


s 


dt 


R 


• 099 


(2 - 49 ) 


In terms of the partial freezing index F 


t 

i-r 


v g dt and realizing 


n-1 n-1 

that R*. = (5 R) + i R , the expression may finally be written as: 

— 2 n 

R 


i* ■» U IL 
n n 1 


* R 2 * 




• • • *(2 « 50) 


For a double-layered system, then, the total freezing index will be: 


F = F x + f 2 


since F_ * XL _1 


F 2 ■ \ ( R i + V 


For a triple-layered system this gives: 


R ( JL + Ug) + R \ ....(2 - 51) 

2 2 2 


E i ( i u i + 


u 2 + V 


+ Vf U 2 ♦ V 


For the general case this beccmes 


+ <§V 

....(2 - 52) 


n _ 

F - > R riu 

k=i k U k 

* t O 

m=k+l m J 

• • • • ( 2 •* 

53) 

This may also be written: 


- 



k-1 



p -r ** + 

E Q 

....(2 - 

5 W 

L ■— * — 

k=l 

m a l — ' 




These last equations are employed in Section 3-02-c in connection 
with the applicability of the nomograph to multilayered systems. In this 
instance ® 1^ <1^ accounting for latent heat effects only, since the 

volumetric heat effects are contained in the effective ^ 


value 


PART III, CONSTRUCTION AND USE OF NOMOGRAPH 

3-01. DESCRIPTION OF THE NOMOGRAPH, - 

A copy of the nomograph as prepared by the authors as part of 
this research is included in an envelope on the inside of the back cover. 
This section of the report -will serve to explain its physical makeup, to 
cite directions for its use, and to discuss its limitations and assump- 
tions, 

The nomograph consists principally of three figures. In addi- 
tion, a table of notation, a list of equations, and an example of its 
use are included on the printed sheet. 

Figure I s Figure 1 is a composite of nomographs already pub- 
lished (20, plate 16 ) for the determination of the volumetric heat of 
the unfrozen and frozen soil C u and Cf 5 and for the evaluation of the 
latent heat of fusion I#. One modification has been made, however, that 
of assuming a value of 0,17 BTU/lb for the specific heat of dry solids 
(instead of 0.20) in the equation for C. The value of 0.17 appears to 
represent a more realistic figure for soil temperatures near the freez- 
ing point (21, page 71 )« 

Figure 2 % Figure 2 presents curves for determining the thermal 
conductivity of the frozen soil, k^, and unfrozen soil, k^. These curves 
are the results of investigations made at the University of Minnesota 
under the immediate supervision of Professor Miles S. Kersten (21, pages 
86, 87, 88, 90) m The writers have changed the form of presentation from 
that originally given but the conductivity values, which are dependent 
primarily on water content and dry density, have not been altered. 

■a) 




Figure 3 g The Depth cf Penetration part of the nomograph may 
he considered to consist of three components? (1) the portion above the 
0( > ydAyy X curves solves the equations for the fusion parameter^^, and 

the thermal ratio 0( 5 (2) the nest of curves represents the relation® 

/ 

ship between and the correction coefficient 5 (3) "the portion 
to the right of the curves solves the equation for the depth of penetra- 
tion X* Origin of the equations and o( s/JL* curves is presented in 
Appendix A* 

Table of Notation s The symbols used on the nomograph are 
generally consistent with prevailing notation,. There ar9 two important 
exceptions., however? (l) *n w has been used for the surface transfer co- 
efficient and (2) n v o n 9 which commonly refers to the mean annual tempera- 
ture., is the number of degrees F by which the mean annual temperature 
exceeds the freezing point of soil moisture® For example., if the freez- 
ing point of soil moisture is taken equal to 3^°^ then,, v 0 35 (mean 
annual temperature) - 31 °<» 

List of Equations ? All equations given in this list are either 
well known or are derived in Appendix A 0 The nomograph is primarily a 
graphical solution of these seven equations® 

Example of the Use of the Frost Penetration Nomograph ? This 
example is discussed in some detail under the following section® 

3-02. DIRECTIONS FOR THE USE OF THE NOMOGRAPH* 

The following paragraphs present concise directions for the 
use of the frost penetration nomograph® These directions are based on 
the writers* present state of knowledge regarding the effect of important 
variables® Limitations and assumptions are given in the following section® 


a* General s - Before the nomograph can he used the following 
terms must first he calculated or estimated. 

(1) Climatic conditions of the locality; v , F, and t* 
These values may he determined directly from weather bureau records or 
estimated from contour maps. (For example. Reference 22, page li+9)« 

(2) Type of Surface; Until further data are available, 

use n ■ 0.9 for all pavement surfaces. The basis for this assumption 

is discussed in Section 3-03-d. 

(3) Soil Properties s For each soil stratum within the 
depth subject to freezing; w, and the soil type (whether predomi- 
nantly sandy or clayey). A soil which is not distinctly either type may 
be treated as explained in Section 3~03*“a. 

b. Homogeneous Soil . - This case is represented by the "Example 
of the Use of the Frost Penetration Nomograph". One can enter, with a 
soft pencil, data for the case at hand in the blanks provided in the table. 
On the nomograph, an arrow pointing into a scale indicates a value which 
is entered. An arrow pointing away from the scale represents a number 
which is read out for use elsewhere. Each line connecting nomographic 
scales, has a sequence number and two arrows. The latter indicate the 
direction of travel required to obtain the answer for that step. 

c. Stratified Soil . - The following steps for this solution 
are suggested: 

Step 1 Determine the approximate depth of frost penetration 
X from Figure of this report. 

Step 2 Determine the approximate correction factor^ from 

Figure 9“A of this report. Use a weighted average 

value of water content within the estimated depth 
of freezing X. 
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A. SEMI-EMPIRICAL CURVES FOR CORRECTION FACTOR X 
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Step 3 Determine L for each stratum from Figure 1 of the 


; nomograph* • ' ’• 

Step I 4 . Determine k for each stratum from Figure 2 of the 


nomograph* 


K 


Step 5 Compute the effective from the formula 


(-) 

k eff 
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+ _ n 


^ 6 A 0 


n 2 _i 

Step 6 Use any effective value for k and determine an ef- 
fective L» Suggestions use k = 1*0 S then 

L = ^eff 0 from Step ^®) 


Step 7 Enter the portion of Figure 3 to the right of the 
0( 9 yiO s X curves on the frost penetration nomo- 
graph with values determined in Steps 2 and 6 to 


obtain a corrected depth of frost penetration X, 

If X differs appreciably from that determined in Step 1, thens 
Step 8 Adjust the value of \ from Step 2* 

Step 9 Correct the value of (— ) from Step 5© 1 1 v ■ 

k eff. 

Step 30 iedetermine X as in Steps 6 and 7® 

Repeat the above steps 3 if necessary* 



d. Negative v 0 (Case where mean annual temperature is below 
the freezing point of soil moisture)* - For the homogeneous case, pro- 
cede with steps 1, 5 » 6* 1 $ and 8 of the "Example of the Use of the Frost 
Penetration Nomograph", then use A® 0*9 to complete the computation 
with steps 13, lii and 15* For the stratified case, useA* 88 0 o 9 ^ n 
Step 2 of Section 3“03-c* 

3-03. LIMITATIONS, ASSUMPTIONS, AND OTHER SPECIAL CASES 0 - It is 
obvious that the accuracy of the nomograph for predicting the depth of 
frost penetration depends to a large extent on the dependability of the 
data used in the computation* This situation is not an unfamiliar one 
in soil engineering problems. If the soil properties of water content 
and dry density are at best only crude approximations or if the climatic 
conditions are not too well known then the nomograph should probably be 
disregarded. In this case, simplified procedures utilizing curves such 
as those in Figures 9-A and 9-B are indicated. 

The simple curve of Figure 9“® which represents a first axjproxima- - 
tion, may be expected to yield predicted depths of frost penetration 
below pavements having granular base courses, differing not mere than 
50 per cent from the actual if the freezing index is greater than 500 
degree days. A second stage approximation may be made by utilizing the 
semi-emperical curves of Figure 9~A to determine a value of the correc- 
tion coefficient A. for use in the equation 

x = 

- 50 - 


0 


■which can be solved by the portion of the nomograph to the right of the 
(X , X curves® The contractors feel that Figures 9-A and 9-B 

might very well be added to the nomograph at some future date* 

The frost penetration nomograph is, of course, subject to all of' 
the limitations and assumptions used to derive the rational formula* 
These are stated and discussed in Section 2-02-e. The importance of 
these assumptions as a whole can best be shown by statistical ^studies 
comparing predicted and actual depths of frost penetration. The reader 


is referred to Section 2-02~d-(2). 

a. Limitation on k* - According to Professor Miles Kersten 


( 21 ) 


the thermal conductivity values from Figure 2 of the nomograph are good 
to + 25 per cent* 

■ If the soil under consideration is a well graded one which is 
neither truly sandy nor a silt or clay then the conductivity may be es- 
timated by interpolation between the dashed and solid curves for the 
given water content and dry density. No specific rules can be offered 
at this time for the interpolation. To quote Professor Kersten, W I my- : - 

self sometimes study the test results on the individual soils in Ap- 

( 21 ) 

pendix 1 to find a soil similar to the one for which k is desired* " 

b. Allowance for Frost Heave * - If freezing temperatures pene- 
trate into a frost susceptible soil which is part of an open system, ice- 


lenses may form and cause considerable heave of the ground surface* The « 
maximum depth of frost penetration in this instance will be smaller than 
for a corresponding case with no ice segregation. In effect, the water 
content of the soil is increased as water flows upward to the zone of 



y 


freezing* This, in turn, increases th9 latent heat of fusion -which is 
the most important soil characteristic affecting the depth of frost peme- 

o 

6 i /l~ 

tration. Since X'vy it follows that the maximum depth of frost pene- 
tration will be smaller* 

If the amount of frost heave due to lease growth can be reason- 
ably estimated, an equivalent average water content, greater than the 
normal, and a corresponding dry density, smaller than the normal, may 
be determined for the soil in the frozen state* These values may then 
be used in the nomograph to predict the depth of frost penetration* 

c* Thaw , - The nomograph at the present time is not directly 
applicable for use in predicting the rate at which thawing occurs during 

i 

the spring months* The writers will study this case under a Corps of 

Engineers contract during the Fiscal Year 1953-i95if* 

( 23 ) 

Carlson and Kersten have shown that the depth of thawing 

below pavements in Alaska may he accurately predicted assuming A, = 1,0 


in the formula 


12 


X\j- 


l*8knl 


where n * l,!*. and I represents the air thawing index* 

d* Surface Transfer Coefficient, n» » The writers have selected 
for use at this time a value of n equal to 0*9 for all pavement surfaces* 
There can be little question, however, that n for bituminous concrete 
pavements is smaller than for Portland cement pavements because of the 
increased absorptivity for solar radiation in the former. The Frost Ef- 
fects Laboratory has, in fact, used values of n equal to 0*75 and 

0*90, respectively, for the two pavement types. These values were 
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arrived' at through an analysis of thermocouple measurements in the two 


of pavements at three site locations. Nevertheless, results of 

.i 

recent statistical studies by the writers and the Frost Effects 
tory using the "modified Berggren formula. Section 2-02-d- 


N V 


that computed depths of freezing are generally too small when pavement 




freezing indices equal to 0,75 and 0,9 are used. Furthermore, there is 
less scattering of statistical results when the air freezing index is 

• f ± • c 

employed. Therefore, it is recommended that until further data and 


analyses are available, a value of n equal to 0.9 be used. 

•-"v:. The surface freezing index for other surface types including- , ; : ,;r 
trees, moss, grass, loam, etc., are undoubtedly considerably smaller 

i 

than 0*9. The writers will not venture a recommendation ’ at this time 
for these Surface types. Some data are available in Reference (l6). 

• A - The following factors are known to affect the surface tempera- • 
ture transfer phenomenons v 

(l) Ambient air temperature; 

: v : (2) Surface connective heat transfer coefficient! 


(3) Solar radiation, direct and diffuse 5 e-y 

(I4.) Absorptivity for Solar radiation; 

(5) Lon g-wa ve radiation (a function of cloudiness and- 

humidity); ' ‘ ' •>•;:. : 

(6) Emissivity and absorbtivity for long-wave radiation 


(7) Precipitation 


, l : j'..' 


The writers- are currently investigating the f actor s listed above- 
and expect to report on their importance as a part of researches for the 


Corps of Engineers during the period 1953-195U® 
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Limitations due to Stratification # - A procedure for apply- 
ing the nomograph to layered soils has been given in Section 3-02-c* Thi 
technique should be restricted to cases which lie between the extremes 
of very small and very large differences in thermal properties from layer 
to layer# In other words, if only very slight variations occur in soil 
type, dry weight and water content, simple average soil constants may 
be determined and the nomograph used as detailed in Section 3-*02-b# On 
the other hand, if the stratification occurs with extreme changes in 
properties, the procedure of Section 3”02-c will lead to only very ap- 
proximate results, which, however, should be as reliable as those com- 
puted by any other technique now available# 

More accurate determinations of this last case would require 
use of the more elaborate methods outlined under Section 2-01-f~(ij.), in 
which event, the nomograph would be of little use# If such procedures 
are employed, the number of depth nodes (the range of n) should be at 
least twice the number of distinct layers# 

3-Olj.# CONCLUSIONS AND RECOMMENDATIONS# - The nomograph as presented 
in the previous sections is felt to be a useful and rapid computational 
aid# 

It is recommended for users of the nomograph that a record be kept 
on the A. - oi, — AO curves of all values of X, obtained from actual 
depth measurements, whenever sufficient data are available for such 
checks# This permits the nomograph to be used in a semi-empirical 
fashion for more reliable predictions in a given locality# 



It is further recommended that a table of n values for different 
surface types be included on the nomograph -when more reliable informa- 
tion becomes available* In addition, it is suggested that curves similar 
to those shown in Figure 9 he- added to the nomograph. These additions 
along with thaw considerations will further the objective of making the 
nomograph a complete entity for determining the depth of freezing and 
thawing* 


PART IV* THAW-CONSOLIDATION PROBLEM 


U-Ol. SYNOPSIS* - 

In the following section* the physical nature of the so-called 


thaw-consolidation problem is discussed* This is followed by a mathema- 


tical formulation of the problem and a discussion of solution techniques. 
It is apparent that a rigorous analytical solution of this 


problem is impossible without questionable simplifying assumptions, 

• r I . . • . - . . ... ....... 

A.1 1,7’ r ;' v \4» v ! v..* /- .. * r ••• _ ' •" ■■■: '• ' . • ** ... 

Numerical methods by hand or I.B.M, and hydraulic models offer the best 


possibilities for practical solutions. Important considerations bearing 
upon the solution of this probleja are discussed in earlier sections, 
principally under Practical Methods of Computation (Section 2-0 1-f) and 


Depth-Time Curves (Secti on 2-02-f ) , 


It is recommended that this problem continue to be explored in 


a limited way with additional I.B.M, and hand numerical solutions but 

p . *. .1 . i . ^ 

that the bulk of the work be postponed until an hydraulic analog is 


constructed along the principles outlined in Section 2-01-f-(l). 

. ;< fv * ’ j * • - ! 

‘ .. \ } , ; • * : ’ . . : n • * : ’ ■ .• 

U-02. PHYSICAL NATURE OF THE PROBLEM. - 


If below-freezing temperatures penetrate into a frost-susceptible 


soil, significant ice segregation is likely to occur. In this state, the 


soil stratum contains an excessive amount of water. More specifically. 


the soil exists at a higher average void ratio and water content than it 


normally would under the weight of the overlying material. 

In the spring, as warm weather approaches, thawing progresses 


from ground surface downward into the supersaturated soil. In effect, 
for each unit of depth the thawing releases a charge of water which was 


stored in the soil in the form of ice lenses® This excess water must , 

- . ’ • \ ’ 

initially flow vertically upward* if a one -dimensional case is visualized* 
through the thawed soil to a drainage surface® The rate at which the 
dissipation occurs depends principally on the permeability of the over- 
lying soil and the distance the water must travel to a drainage surface® 

After the soil stratum has completely thawed, the excess water is eventu- 
ally entirely dissipated through a seepage diffusion process® The soil 
has adjusted, then to a stable void ratio under the weight of the overlying soil. 

J+-03® MATHEMATICAL FORMULATION ® - 


Before formulating the equations and boundary conditions applying 
to this problem it is important to investigate the pressure-void ratio ' 
relationship of an element of soil during the thawing process® 

In the classical one-dimensional consolidation theory a straight 
line relationship between void ratio and pressure is assumed® It is 
apparent that for frost-heaves of appreciable magnitude the void ratio of 
the frozen soil will increase to a value far greater than that given by 


the usual straight line relationship extrapolated to zero pressure® The 

straight line relationship thus needs clarification® 

let e and pp> Figure 10, represent the void ratio and inter- 

granular pressure, respectively, in an element of soil before ice-segregation 
occurs® Point A thus represents an equilibrium state under the weight of 
the overlying soil® Line CD through point A represents a portion of the 
pressure-void ratio relationship of the element® Its local slope*^?. » is 
thus the coefficient of compressibility a^* 

Hhen significant ice segregation occurs, the void ratio in effect 
increases to some value e..® Most of the load is actually carried by the ioe. 







I. suddenly thawed, the intergranular pressure would thus be very nearly 
zero as the element begins to consolidate* Assume, then, that point B 
represents the soil both in the frozen state and immediately after thaw* 
After thawing, the element follows some path from B to A as it 
consolidates under the weight of the overlying soil* The path is probably 
given by a curve similar to that indicated in Figure 10 by the heavy dashed 


line * 


As a first approximation, it is reasonable to assume that the 


element follows the path BC'A, Under this assumption the physical seepage 

and consolidation process may be described as follows. On first thawing 

the element has a void ratio e^, and a hydrostatic excess water pressure 

equal to p^« The excess water pressure remains constant in the element 

until enough stored water has bean dissipated to lower the void ratio to 

Flow occurs throughout this period under a constant excess water 

pressure p^ since the intergranular pressure is zero. In consolidating 

from e 7 to e 0 the soil follows C ? A (with a = constant, determined at A) 
o d v 

as the excess pressure is dissipated according to the unsteady flow 


I ■* Cr V 


diffusion process, governed by the fundamental equations 

= k 4.5 
d x , . 

i| + « 0 

$t d x 

As = — A 


o 9 




1 f. 




0000 


• o « • 


1+e 

Au = K &h = A 


(U~2) 


(k™3) 




e 


• « O' 


.(U-W, 




which may be combined into the classical diffusion equation for the 

S-b fi’cc i crif : . ;i. \ 2 ... V • • * , - - .-ss £ ' 

2 a *3— ii - , , 0 . (U"5) 

O * d t 


J.’. ‘ 


hydrostatic excess pressure u, as. 

k (l+e) 


where c- v ” 


a v 
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U-Olu SOLUTION TECHNIQUES. - 

; - • 

Even with the simplifying assumptions as stated above, the thaw- 
consolidation problem as formulated in the previous section is not amenable 
tc closed form analytical solution. 

This difficulty suggests the application of methods similar to 
those outlined in Section 2-01-f these are detailed in the sections to follow, 

a. Hydraulic Model . - The use of an hydraulic model represents 
of the most efficient methods of solving the thaw-consolidation problem. In 
addition, the model presents a graphic picture of the state of consolidation 
within the soil stratum at any time. The relationship of an hydraulic model 
to the thermal diffusion problem with latent heat is described in Section 
2-01-f-(l). Its relationship to the consolidation process is presented in 
References 3 and Only the direct application of the model to the thaw- 
consolidation process will be discussed here. 

A simplified hydraulic model for the one dimensional thaw- 
consolidation process is shown in Figure 11. This particular model is 
based on the assumption that the soil is homogeneous and that each soil 
element follows a path similar to line BC*A, Figure 10, as it. consolidates 
following thaw. It is further assumed for this discussion that the co- 
efficient of permeability k remains constant during the consolidation process. 

As a result of ice-segregation in the frozen portion of 
the soil an excess amount of water is stored in each layer. This excess 
water is represented by the volume of water in the reservoir and corre- 
sponding standpipe. The water volume in the reservoir represents the void 
ratio change e^ - e^. Figure 10, while that in the standpipe corresponds to 



in each soil element 


e 3 " © 2 * After thaw, the initial excess pressure 
is equal to the effective weight of the overlying soil* This assumption 

fixes the initial water level in the standpipes* 

■When thawing penetrates a distance A x, the excess water 
stored in the first layer of soil is released* In essence, valve A is 
opened and steady flow commences toward the adjacent tank I under a head 
differential equal to u^* As soon as the reservoir is empty unsteady 
flow from the standpipe occurs under continously decreasing head* 

Successive valves are opened as each layer thaws* When the 
depth of thaw equals the soil within the unfrozen lower portion of the 
stratum takes part in the diffusion process and swells as water flows ver- 
tically downward toward the bottom drainage surface* On the model this 
condition is initiated by opening the last valve E* 

Any specific thavf-consolidation problem may be treated 
using two separate hydraulic models* The first model is set up for the 
thermal problem of thaw while the second, similar to that in Figure 11, 
represents the accompanying consolidation process* The first model in- 
dicates the time when each valve in the second should be opened* 

It is theoretically possible to construct the model such 
that it can represent complex pressure-void ratio functions during thaw, 

variations in permeability with void ratio, and other features of stratifi- 

! 

cation and non-homogeneity* These refinements may be incorporated when 
knowledge regarding the actual thaw-consolidation characteristics of soil 
has been advanced* Laboratory and full -scale field tests are needed for 
this purpose* 
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b. Numerical Methods. - Numerical methods of solution 


applicable to the thaw-consolidation problem parallel the methods dis- 
cussed in Sections 2-0 1-d, 2-01-f-(3) and 2-01-f-(U) . They can be 
divided into two -types, hand computations and machine computations* 

(1) Hand Computations. - The technique of solving a 
problem of this type is best illustrated by an example. Such an exampl# 
is presented in Appendix C* 

(2) Machine Computations. - The use of automatic 
digital computing equipment such as the I.B.M. Card Programmed Calculator 
represents a rapid but relatively costly method of obtaining solutions 

to the thaw-consolidation problem. Techniques very similar to those 
described in Section 2-01-f-(3) would be used* 
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Appendix A 


DERIVATION OF A RATIONAL FORMULA FOR THE 

Prediction of frost penetration 


A-Ol. BASIC ASSUMPTIONS. - 

We seek a solution to the thermal problem of a semi -infinite soil 
mass of uniform properties and at a uniform initial temperature subjected to 
changes of temperature at the surface, which are assumed uniform over the 
surface extent to yield a one-dimensional problem. 

Figure h illustrates the nomenclature and significant variables 
in this situation. Farther assumptions are best listed in the form of 
specific conditions as follows: 

Condition I - AT THE GROUND SURFACE 

It is assumed that the surface temperature is suddenly changed 
from an initial temperature v Q above freezing to a temperature r g below 
freezing. This temperature value is then maintained constant and uniform 
over the entire surface. 

Condition II - IN THE FROZEN SOIL 

It is assumed that the diffusion equation: 


a f V v 


2 _ _ d 

f 


d t 


with a^=k^/c^, measuring the diffnsivity of the frozen soil is satis- 
fied throughout the frozen soil mass, subject to the surface temperature 
condition (I) and the latent heat condition (ill) at its boundaries. 
Condition III - AT THE MOVING FROST INTERFACE 


It is assumed that at the interface between the frozen soil 
(above) and the unfrozen soil (below) the temperature remains constant 


at the freezing point of the soil moisture. It is further assumed that 
the heat flow upward just above the interface in the frozen soil is equal 
to the sum of the heat flow just below the interface in the unfrozen soil 
plus the heat flow due to the removal of the latent heat of fusion of the 
soil moisture as it freezes. 

Condition IV - IN THE UNFROZEN SOIL 


It is assumed that the diffusion equations 


1 V 

u y 




U 


U 


d t 


with a,^ = measuring the diffusivity of the unfrozen soil, is sat- 

isfied throughout the unfrozen region, subject to the latent heat condi- 
tion (III) at the frost interface (the upper boundary of the unfrozen 
soil) and to the lower boundary condition that, at all times, the tempera- 
ture approaches the initial temperature for sufficiently great depths. 


A-02. GENERAL SOLUTION. - 

Following the notation of Figure I 4.9 the solution for v^, which 
satisfies conditions I and II may be written: 


- V S + A arf (— ) 

2 ^ 


O • 0 (A— 1 ) 


while the solution for v u satisfying condition IV may be written: 


v *= v + B 
u o 


1 - 


erf (. 


•) 


• c • o ( ) 


2 W 

At the interface (x ** X), in order to satisfy condition III, 

where the freezing point is taken arbitrarily as zero degrees, it is nec- 
essary that 


v f = v u 


•«•« (A— 3 ) 






\ 


A-2 


from which: 


A erf ( 


B 


[*- 




= v 


s 


erf ( 


2\/vE 


) “ - t , 


. . • • (A-I4.) 


• • • • (A-5 ) 


But since these last conditions must be satisfied for all values of time. 


then 


or 


X / = constant 


X\ft 


..•.(A-6) 


• • • • (A-7 ) 


Thus the constant depends upon v g , v^ and the thermal coefficients* 

The manner of this dependence may be found by considering the latent heat 
requirement of condition III. This thermal continuity condition relates 
the temperature gradients each side of the interface to the rate of move- 


ment of the interface in the form: 

... 


*11 


dX 


d* 


dt 


... • (A-8) 


which becomes, upon substituting the expressions (A-I4.), (A-5) and (A-7), 


performing appropriate differentiation, and simplifying: 


Vf e 


-t 2 /^ 


v k 
o u 


e 


- 


, L 


'fruTf. ' * p - erf (f/2^) 




2 


• • • • (A-9 ) 


Furthermore, by making the substitutions 

T 0 V\ c t 

$ = V k f c iAu c f 


....(A-10) 


....(A-n) 



su- = V s c f/ L 


...(A-12) 


Z = K/Z 


o o o o 


(A-13 ) 


it is possible to rewrite Equation (A° 9 ) in a simplified non-dimensional 


form, namely 




e 




J> 2 Z 2 


_erfZ S (l - erfgZ) J 


= V /Tf 


...(A-lU) 


If all the terms in Z are then carried to the right hand side. 


one has a direct relation between/^ and C>f , § , and Z in the form 

^ -b z 


- z 2 


>*/<K5,n -y/Tf 2 / 


erf Z 


(1 - erf 5 Z )_ 


»...(A-15) 


This last nay then be inverted to obtain a transcendental relation for 2 
in terms of the parametersQf', and Jjl 5 which may be solved graphically 


to give, in symbolic form 


Z “ Z W.Si/) 

Thus using Equation (A-13), one may obtain 
X s 2 ^"a^ • Z 


0 6 0 0 


(A- 16 ) 


OOOO 


(A- 17 ) 


and, finally, for the frost penetration depths 

X = 2 Z Ja f t 000 . (A-18) 

It is found convenient to rearrange this last expression so that the 


radical is expressed in terms of^CL in the forms 


2 Z 

^u- 




OOOO 


(A-19) 



By defining the new d imensi onless correction coefficient X as 

A.=\/tT- = ....(A-20) 


and by substituting under the radical in Equation (A- 19 ) for the physical 


variables, one finally obtains for X, the general expression 




••••(A-21) 


for time measured in days # The correction coefficient is, again, a 
function of the three dimensionless parameters^, £ » yOL • 

A-03. PHYSICAL SIGNIFICANCE OF TERMS AND PARAMETERS. - 

a. Surface Temperature Change (v ). - In order to apply Equa- 

g 

tion (A-21) to act-al cases in which the surface temperature does not 

change suddenly as assumed in this derivation, the average temperature 

below freezing during the freezing period is taken as v , Thus, in terms 

s 

of the freezing index F (in degree-days) and the freezing period t (in 


days), the surface temperature v„ may be written 


v = F/t 
s 


• • • • (A-22) 


or, solving for F, 


F = v t 
s 


...(A-23) 


The effects of the actual shape of the freezing index vs» time 
curve — as compared with the assumed linear shape — are not significant 
so long as the freezing period is not too short nor the freezing index too 
small (i.e., cases in which the depth of freeze would normally be fairly 
small.) This problem is treated in Section 2-02-f. 

b. Thermal Ratio (&( ), - The thermal ratioO^, defined as 

v 0 C tAs C f ••••(A-2U) 


A-5 



measures the ratio of heat stored initially in the unfrozen soil to the 
boat loss in the frozen soil c If it may be assumed, as with the existing 
formulas, that the difference between the volumetric heats and is 
not usually significant, the thermal ratio^may be written 

& = v oA s ••••(A-25) 

which is the ratio of th9 initial ground temperature (or mean annual tem- 
perature) above the freezing point to the average surface temperature be- 
low freezing during the freezing period. 

In terms of the freezing index F and the freezing period t, 
this may then be written 

OC « v 0 t/F ... . (A~26) 

Co Pi f f usivity Ratio ( 5 ) . - The (root) diffusivity ratioS, 

defined a.= 

^ = V a f/ a u o««e(A-27) 

measures th9 relative values of diffusivity a = |?/C in th6 frozen and 

unlV^^r: soils. It is clear that for most soils of low moisture content 
S is approximately unity. Tabulated herewith are typical values of § for 
representative soil types and moisture contents. 


TABLE A - I 


TYPICAL VALUES OF THE DIFFUSIVITY RATIO 
Assumed Dry Density * 110 Ibs/ft^ 


Moisture Content 

% 


Values of 5 



Silt or Clay 

Sand 

0 

1.00 

1.00 

5 

1.07 

0.92 

10 

1.15 

1.14 

15 

1.22 

1.31 

20 

1.50 

1.50 


Source of Data for k and C: "Final Report, Laboratory Research for the 

Determination of the Thermal Properties of Soils" C, of E«, St. Paul 
District, June 19^*9 • 


The effect of variations in S may be found directly from the 
A-<x - /CL graph. Figure 5 > by making use of an equivalent 0( -value given 


by the expression: 




c<_ 

5 


1 - erf, 
_ 1 - erf 




-J 


- (S 2 -i) X 2 ^ 


, (A-28) 


•••• (A-29 ) 


erf5X\^¥ 

which can be approximated by the empirical equation: 

5 ~ 5 L 2 

A comparison of exact values of determined from Equations (A-15) and 
(A-20) with estimated values using Figure 5 and Equation (A-29) shows only 
negligible differences for all practical values of § . 


A-7 


The following table indicates the relationship between 0 ( , 5 
and as given by Equation (A-29)« 


TABLE A - II 


EQUIVALENT THERMAL RATIO <*fc> 
Values of 4c - £ST fl + 1 

^ 5 L 




VALUES OF 

<x 

VALUES OF § 

0.9 

1.0 

1.1 

1.2 

1.3 

1.U 

1.5 

0 

0 

0 

0 

0 

0 

0 

0 

l 

1.07 

1.00 

0.9k 

0.89 

0.85 

0.82 

0.79 

2 

2.16 

2.00 

1.87 

1.76 

1.67 

1.59 

1.53 

3 

3*25 

5.00 

2.80 

2.62 

S.I 4 S 

2.36 

2.25 

k 

U.35 

It. 00 

3.72 

3.47 

3.28 

3.11 

2.97 


Since it has been demonstrated that a 1% error in water content 
produces only about a 1.5$ change in the value of , the above tabulation 
in conjunction with Figure 5 would indicate less than a 1 $ change in the 
computed depth of oenetration due to this variation,, Moreover, typical 
soils have § values of the order of 1.15, which would indicate an effect- 
ive value of which is roughly 15 % smaller than that computed as- 

suming § = 1*0, yet this effect can produce a variation inland there- 
fore a change in predicted depth of less than 5$ for typical thermal con- 
ditions. 

Since inclusion of a tern involving S would increase the com- 
plexity of representation of the formula in graphical form, and would 
further increase the tendency to predict depths of penetration which are 













too deep, it seems reasonable to base at least preliminary calculations 
on the assumption that S B 1*0* 


Consistent, then, with the previous assumption that C 


C u ■ C f 


and 5 83 /(—£—) (_1L_) = 1.0, there would follow then the necessary condi - 

* *Si ®f 

tion that k ° k^ ■ k^ • In practical applications, this would suggest 
the use of average values of C and k. 

d. Rision Parameter (/Q . - The fusion parameter/^ , defined 


as 




s 


C/L 


....(A-30) 


measures the heat removed in the frozen soil (below the freezing point) 
compared to the latent heat of the soil moisture. Ihis characteristic 


may be written in terms of F, t, and C in the form 

= rc/Lt 


....(A-31) 


When//. 38 0, the only significant soil properties affecting the 
depth of frost penetration are the latent heat L, and the thermal conduc- 
tivity, k^j on the other hand* as^/c becomes large, the stored heat in the 
soil volume becomes proportionately more significant. 

e. Correction Coefficient (X ). - If Equation (A-22) is sub- 
stituted into Equation (A-21) we obtain the more familiar expression 



....(A-32) 


and^ , therefore, becomes a correction coefficient applied to the calcu- 
lated depth of penetration due to latent heat only. 

The value ofX ma y be found from Figure 5 which shows the rela- 
tionship between 0( ty /jL and X. for the case of 6 " 1. Data for this figure 
were obtained by assuming values of z and solving Equation (A-15) for jU. . 
One QyU. was obtained^, followed directly from Equation (A-20). 


A-Oiu COMPARISON WITH BERGGREN FORMULA. - 

(19) 

As presented in Section 2-02-b. Shannon has written the 
( 18) 

Berggren formula ; in the form: 

X = 2B \J a f t ....(A-33) 

By comparison with Equation (A-18) of Section A-02, it is clear that the 
term B in Equation (A-33) is identical with the term Z in Equation (A-18)* 
Thus the only significant difference between the older Berggren formula 
and the present rational formula lies in the inversion of the nuclear 
radical term to permit expressing the correction factor^. to the depth 
of frost penetration that would occur if only latent heat were present, 
rather than the correction term B — Z applied to the depth occurring if 
only volumetric heat is significant* Since for all practical problems 
the effects of L are large compared to those of C (i.e., the yCC values 
are low), this seems to be the more useful form for the same fundamental 
solution. 

For the above reasons, it would seem fitting to refer to the 
rational formula proposed here as the "modified Berggren formula"* 




A-iO 


APPENDIX 3 


TABLES OF I.B.M, SOLUTIONS 


B-Ol. GENERAL. - 

As outlined in Section 2-01-f-(3). p&rt of this investiga- 
tion -tvas concerned with obtaining precise numerical solutions to the 
difference equations which represent the thermal transients in the 
ground during the freezing and thawing processes. These equations are 
presented and explained in the section cited. 

Solutions to these equations were run on the I.B.M. Card Pro- 
grammed Calculator of the M. I. T. Statistical Services Division, v/ith 
values and conditions indicated in the following tabulation. 


B-02. INDEX TO SOLUTIONS. - 


TABLE 

SURFACE 
TEMPERA TUT- 
VARIATION 

PARAMETERS 

RANGE 



COEFF. 

£ 



n 

K 

B-I 

Step 

0.5 

0+* 

8 

2k 

0.25 

B~II 

Step 

0.5 

1.0 

8 

24 

0.25 

B-III 

Step 

0.5 

2.C 

3 

2 h 

0.2^ 

B-rv 

Step 

0.5 

5.0 

8 

2k 

0.25 

B~V 

Step 

0.5 

U.0 

O 

O) 

0.25 

B-VI 

Step 

0.5 

5.0 

n 

2k 

C.25 

B-VII 

Step 

1.0 

C+* 

n 

O 

24 

0.25 

B-VI I I 

Step 

1.0 

1.0 

8 

2k 

0.25 

B-IX 

Step 

1.0 

2.0 

8 

2k 

0.25 

B-X 

Step 

1.0 

3.0 

5 

2k 

0.25 

B-XI 

Step 

1.0 

l*.o 

0 

24 

c.25 

B-XII 

Step 

1.0 

5.0 

8 

2k 

c.25 

B-XIII 

Cyclic 

1.0 

0 

l 

108 

0.1*0 

B-xrv 

Cyclic 

loO 

1.0 

7 

103 

0.1*0 


♦Temperature initially at the freezing point, soil unfrozen. 


B-03. TABULATED VALUES. 


The tabular entries in the tables give the value of e . cor- 
responding to each value of n and k. These values for n = 0 correspond 


B-l 











to the surface temperature variation, while the remaining values repre- 
sent subsurface temperatures. The transformation of these dimension- 
less values, together with the time measure (k) and depth measure (n), 
are explained in Section 2-01-f-(3)» 


B-2 


TABLE B-I 

>u * 0.5 <* = 0+ 


k 

n VALUES 

VALUES 

0 

1 

2 

3 

4 

5 

6 

7 

8 

0 

-0.50 

0 

0 

0 

0 

0 

0 

0 

0 

1 

-o.5o 

0 

0 

0 

0 

0 

0 

0 

0 

2 

-o.5o 

0 

0 

0 

0 

0 

0 

0 

0 

3 

-o.5o 

0 

0 

0 

0 

0 

0 

0 

0 

4 

-0.50 

0 

0 

0 

0 

0 

0 

0 

0 

5 

-o.5o 

0 

0 

0 

0 

0 

0 

0 

0 

6 

-0.50 

0 

0 

0 

0 

0 

0 

0 

0 

7 

-0.50 

0 

0 

0 

0 

0 

0 

0 

0 

8 

-0.50 

0 

0 

0 

0 

0 

0 

0 

0 

9 

-0.50 

-0.12 

0 

0 

0 

0 

0 

0 

0 

10 

-o.5o 

-0.19 

0 

0 

0 

0 

0 

0 

0 

11 

-o.5o 

-0.22 

0 

0 

0 

0 

0 

0 

0 

12 

-0.50 

-0.23 

0 

0 

0 

0 

0 

0 

0 

13 

-OoO 

-0.2U 

0 

0 

0 

0 

0 

0 

0 

1 h 

-0.50 

-0.25 

0 

0 

0 

0 

0 

0 

0 

15 

-0.50 

-0.25 

0 

0 

0 

0 

0 

0 

• 0 

16 

-0.50 

-0.25 

0 

0 

0 

0 

0 

0 

0 

17 

-0.50 

-0.25 

0 

0 

0 

0 

0 

0 

0 

18 

-0.50 

-0.25 

0 

0 

0 

0 

0 

0 

0 

19 

-0.50 

-0.25 

0 

0 

0 

0 

0 

0 

0 

20 

-o.5o 

-0.25 

0 

0 

0 

0 

0 

0 

0 

21 

-o.5o 

-0.25 

0 

0 

0 

0 

0 

0 

0 

22 

-0.50 

-0.25 

0 

0 

0 

0 

0 

0 

0 

23 

-0.50 

-0.25 

0 

0 

0 

0 

0 

0 

0 

24 

-0.50 

-0.25 

0 

0 

0 

0 

0 

0 

0 


TABLE B-l 











TABLE B-II 

yU = 0.5 «* 5 1.0 


k 

VALUES 

n VALUES 

0 

1 

2 

3 

4 

5 

6 

7 

8 

0 

-0.50 

0.50 

0.50 

0.50 

o.5o 

0.50 

o.5o 

0.50 

o.5o 

1 

-0.50 

0.25 

o.5o 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

2 

-0.50 

0.12 

o.l*l* 

0.50 

0.50 

o.5o 

0.50 

0.50 

0.50 

3 

-0.50 

0.05 

0.38 

o.48 

0.50 

0.50 

0.50 

0.50 

0.50 


-o.5o 

0 

0.32 

0.48 

0.50 

o.5o 

o.5o 

0.50 

0.50 

5 

-0.50 

0 

0.28 

0. U6 

0.49 

0.50 

0.50 

0.50 

o.5o 

6 

-0.50 

0 

0.26 

0.1*6 

0.49 

0.50 

0.50 

0.50 

0.50 

7 

-0.50 

0 

0.2lt 

0.46 

0.1*8 

o.4o 

0.50 

0.50 

0.50 

8 

-0.50 

0 

0.24 

0.45 

0.1*8 

0.49 

0.50 

0.50 

o.5o 

9 

-0.50 

0 

0.23 

0.45 

0.48 

0.49 

0.50 

0.50 

0.50 

10 

-o.5o 

0 

0.23 

0.45 

0.47 

0.49 

0.50 

0.50 

0.50 

11 

-o.5o 

0 

0.23 

0.45 

0.47 

0.49 

0.50 

0.50 

o.5o 

12 

-o.5o 

0 

0.22 

0.45 

0.47 

0.48 

0.49 

o.5o 

0.50 

13 

-0.50 

0 

0.22 

0.1*1* 

0.47 

0.48 

o.49 

o.5o 

0.50 

14 

-0.50 

0 

0.22 

0.1*1* 

0.46 

0.48 

0.49 

0.50 

0.50 

15 

-o.5o 

0 

0.22 

0 

0.1*6 

0.48 

0.49 

o.5o 

0.50 

16 

-o.5o 

0 

0.22 

o.l*l* 

0.46 

0.48 

o.49 

0.50 

o.5o 

17 

-0.50 

0 

0.22 

o.l*l* 

0.1*6 

0.48 

o.49 

0.49 

0.50 

16 

-0.50 

0 

0.22 

o.l*l* 

0.1*6 

0.48 

0.49 

0.49 

0.50 

19 

-0.50 

0 

0.22 

o.l*l* 

0.1*6 

0.47 

0.48 

0.49 

0.50 

20 

-0.50 

-.06 

0.22 

0.104 

0.1*6 

o.47 

0.48 

o.49 

0.50 

21 

-0.50 

-.10 

0.20 

o.lO* 

0.1*6 

0.47 

o.48 

o.49 

0.50 

22 

-o.5o 

-.12 

0.19 

o.43 

0.46 

0.47 

0.48 

0.49 

0.50 

23 

-o.50 

-.11* 

0.17 

0.43 

o.45 

o.47 

0.48 

0.49 

0.50 

24 

-o.5o 

-.15 

0.16 

0.42 

0.45 

o.47 

0.1*8 

0.49 

0.50 


TABLE B-II 










table B-III 

/jl = 0.5 00 ~ 2*0 



TAE1£ B-III 









TABLE B-IV 
* 0.5 06 * 3.0 


n VALUES 

0 

#H 

2 

3 

1 + 

5 

6 

7 

8 

-0.50 

1.50 

1.50 

1.50 

1.50 

1.50 

1.50 

1.50 

1.50 

-0.50 

1.00 

1.50 

1.50 

1.50 

1.50 

1.50 

1.50 

1.50 

-0.50 

0.75 

1.37 

1.50 

1.50 

1.50 

1.50 

1.50 

1.50 

-0.50 

0.59 

1.25 

1.1*6 

1.50 

1.50 

1.50 

1.50 

1.50 

-0.50 

0.1*8 

1.1 k 

1.1*5 

1.1*9 

1.50 

1.50 

1.50 

1.50 

-0.50 

o.uo 

1.05 

1.1*3 

1.1*8 

1.1*9 

1.50 

1.50 

1.50 

-0.50 

0.31* 

0.98 

1.1*1 

1.1*7 

1.1*9 

1.1*9 

1.50 

1.50 

-0.50 

0.29 

0.93 

1.1*0 

1.1*6 

1.1*9 

1.1*9 

1.1*9 

1.50 

-0.50 

0.25 

0.89 

1.39 

1.1*5 

1.1*6 

1.1*9 

1.1*9 

1.50 

-0.50 

0.22 

0.86 

1.38 

1.1*1* 

1.1*8 

1.1*9 

1.1*9 

1.50 

-0.50 

0.20 

0.83 

1.37 

1.UJU 

1.1*7 

1.1*9 

1.1*9 

1.50 

-0.50 

0.18 

0.81 

1.37 

1.1*3 

1.1*7 

1.1*9 

1.1*9 

1.50 

-0.50 

0.17 

0.79 

1.36 

1.1*2 

1.1*6 

1.1*8 

1.1*9 

1.50 

-0.50 

0.16 

0.78 

1.36 

1.1*2 

1.1:6 

1.1*8 

1.1*9 

1.50 

-0.50 

0.15 

0.77 

1.35 

1.1*1 

1.1*5 

1.1*8 

1.1*9 

1.50 

-0.50 

O.lli 

0.76 

1.35 

1.1*1 

1.1*5 

1.1*7 

1.1*9 

1.50 

-0.50 

0.1U 

0.75 

1.31* 

1.1*0 

1.1*1* 

1.1*7 

1.1*8 

1.50 

-0.50 

0.13 

0.75 

1.31* 

1.1*0 

1.1*1* 

1.1*7 

1.1*8 

1.50 

-0.50 

0.13 

0.71* 

1.33 

1.39 

1.1*1* 

1.1*6 

1.1*8 

1.50 

-0.50 

0.12 

0.71* 

1.33 

1.39 

1.1*3 

1.1*6 

1.1*6 

1.50 

-0.50 

0.12 

0.73 

1.33 

1.39 

1.1*3 

1.1*6 

1.1*8 

1.50 

-0.50 

0.12 

0.73 

1.32 

1.38 

1.1*3 

1.1*6 

1.1*8 

1.S0 

-0.50 

0.12 

0.73 

1.32 

1.38 

1.1*2 

1.1*5 

1.1*7 

1.50 

-0.50 

0.12 

0.72 

1.32 

1.37 

1.1*2 

1.1*5 

1.1*7 

1.50 

-0.50 

0.11 

0.72 

1.32 

1.37 

1.1*2 

1.1*5 

1.1*7 

1.50 


TABLE B-IV 










TABLE B-V 

/X = 0.5 06 s h.O 


k 

VALUES 

n VALUBS 

0 

1 

2 

3 

4 

5 

6 

7 

8 

0 

-0.50 

2.00 

2.00 

2.00 

2.00 

2.00 

2.00 

2.00 

2.00 

1 

-0.50 

1.37 

2.00 

2.00 

2.00 

2.00 

2.00 

2.00 

2.00 

2 

-o.5o 

1.06 

1.84 

2.00 

2.00 

2.00 

2.00 

2.00 

2.00 

3 

-0.50 

0.87 

1.68 

1.96 

2.00 

2.00 

2.00 

2.00 

2.00 

4 

-0.50 

0.73 

1.55 

1.94 

1.99 

2.00 

2.00 

2.00 

2.00 

5 

-0.50 

0.63 

l.IOi 

1.91 

1.98 

1.99 

2.00 

2.00 

2.00 

6 

-0.50 

0.55 

1.35 

1.89 

1.96 

1.99 

1.99 

2.00 

2.00 

7 

-0.50 

0.4? 

1.29 

1.88 

1.95 

1.98 

1.99 

1.99 

2.00 

8 

-0.50 

0.44 

1.23 

1.86 

1.94 

1.98 

1.99 

1.99 

2.00 

9 

-o.5o 

o.4o 

1.19 

1.85 

1.93 

1.97 

1.99 

1.99 

2.00 

10 

-0.50 

0.38 

1.16 

1.84 

1.92 

1.97 

1.99 

1.99 

2.00 

11 

-o.5o 

0.35 

1.13 

1.83 

1.91 

1.96 

1.98 

1.99 

2.00 

12 

-0.50 

0.34 

1.11 

1.83 

1.91 

1.95 

1.98 

1.99 

2.00 

13 

-0.50 

0.32 

1.10 

1.82 

1.90 

1.95 

1.98 

1.99 

2.00 

14 

-o.5o 

0.31 

1.08 

1.81 

1.89 

1.94 

1.97 

1.99 

2.00 

15 

-0.50 

0.30 

1.07 

1.81 

1.88 

1.94 

1.97 

1.98 

2.00 

16 

-0.50 

0.30 

1.06 

1.80 

1.88 

1.93 

1.96 

1.98 

2.00 

17 

-o.5o 

0.29 

1.05 

1.80 

1.87 

1.93 

1.96 

1.98 

2.00 

18 

-0.50 

0.28 

1.05 

1.79 

1.87 

1.92 

1.96 

1.98 

2.00 

19 

-0.50 

Q.28 

1.04 

1.79 

1.86 

1.92 

1.95 

1.98 

2.00 

20 

-0.50 

0.28 

1.04 

1.78 

1.86 

1.91 

1.95 

1.97 

2.00 

21 

-0.50 

0.27 

1.03 

1.78 

1.85 

1.91 

1.95 

1.97 

2.00 

22 

-0.50 

0.27 

1.03 

1.78 

1.85 

*1.90 

1.94 

1.97 

2.00 

23 

-0.50 

0.27 

1.03 

1.77 

1.84 

1.90 

1.94 

1.97 

2.00 

24 

-o.5o 

0.27 

1.02 

1.77 

1.84 

1.90 

1.94 

1.96 

2.00 


TABLE B-V 










= 



k 

VALUES 


8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 


TABLE B-VI 

0.5 oc= 5.0 


n VALUES 


3 




-0.50 

2.50 

2.50 

2. 5'0 

2.50 

2.50 

2.50 

2.50 

2.50 

-o.5o 

1.75 

2.50 

2.50 

2.50 

2.50 

2.50 

2.50 

2.50 

-o.5o 

1.37 

2.31 

2.50 

2.50 

2.50 

2.50 

2.50 

2.50 

-o.5o 

l.U* 

2.12 

2.U5 

2.50 

2.50 

2.50 

2.50 

2.50 

-0.50 

0.98 

1.96 

2.1*2 

2.1*8 

2.50 

2.50 

2.50 

2.50 

-0.50 

0.35 

1.83 

2.39 

2.1*7 

2.1*9 

2.50 

2.50 

2.50 

-o.5o 

0.76 

1.72 

2.37 

2.1*6 

2.1*9 

2.1*9 

2.50 

2.50 

-0.50 

0.69 

1.61* 

2.35 

2.1*U 

2.1*8 

2.1*9 

2.1*9 

2.50 

-o.5o 

0.63 

1.58 

2. 3U 

2.1*3 

2.1*7 

2.1*9 

2.1*9 

2.50 

-0.50 

0.59 

1.53 

2.32 

2.1*2 

2.1*7 

2.1*9 

2.1*9 

2.50 

-o.5o 

0.55 

1.1*9 

2.31 

2.1*1 

2.1*6 

2.1*8 

2.1*9 

2.50 

-0.50 

0.22 

1.1*6 

2.30 

2.1*0 

2.1*5 

2.1*8 

2.1*9 

2.50 

-0.50 

0.50 

1.1*1* 

2.29 

2.39 

2.1*5 

2.1*8 

2.1*9 

2.50 

-0.50 

0.1*9 

1.1*2 

2.29 

2.38 

2.1*1* 

2.1*7 

2.1*9 

2.50 

-0.50 

0.1*7 

1.1*0 

2.28 

2.37 

2.1*3 

2.1*7 

2.1*8 

2.50 

-0.50 

0.1*6 

1.39 

2.27 

2.36 

2.1*3 

2.1*6 

2.1*8 

• 

2.50 

-0.50 

0.1*5 

1.38 

2.26 

2.36 

2.1*2 

2.1*6 

2.1*8 

2.50 

-0.50 

0.1*5 

1.37 

2.26 

2.35 

2.1*1 

2.1*5 

2.1*8 

2.50 

-0.50 

0.1*1* 

1.36 

2.25 

2.31* 

2.1*1 

2.1*5 

2.1*7 

2.50 

-0*50 

0.1*1* 

1.35 

2.25 

2.31* 

2.1*0 

2.1*1* 

2.1*7 

2.50 

-o.5o 

0.1*3 

1.35 

2.21* 

2.33 

2.1*0 

2.1*1* 

2.1*7 

2.50 

-o.5o 

0.1*3 

1.31* 

2.21* 

2.32 

2.39 

2.1*1* 

2.1*7 

2.50 

-0.50 

0.1*2 

1.31* 

2.23 

2.32 

2.39 

2.1*3 

2.1*6 

2.50 

-o.5o 

0.1*2 

1.33 

2.23 

2.31 

2.38 

2.1*3 

2.1*6 

2.50 

-o.5o 

0.1*2 

1.33 

2.23 

2.31 

2.38 

2.1*2 

2.1*6 

2.50 


TABLE B-VI 









TABLE B-VII 
ya * 1.0 oc = o+ 


k 

VALUES 

n VALUBS 

o 

1 

2 

3 

b 

5 

6 

7 

8 

0 

- 1.00 

0 

0 

0 

0 

0 

0 

0 

0 

1 

- 1.00 

0 

0 

0 

0 

0 

0 

0 

0 

2 

- 1.00 

0 

0 

0 

0 

0 

0 

0 

0 

3 

- 1.00 

0 

0 

0 

0 

0 

0 

0 

0 

k 

- 1.00 

0 

0 

0 

0 

0 

0 

0 

0 

5 

- 1.00 

- 0.25 

0 

0 

0 

0 

0 

0 

0 

6 

- 1.00 

- 0.38 

0 

0 

0 

0 

0 

0 

0 

7 

- 1.00 

-o.U* 

0 

0 

0 

0 

0 

0 

0 

8 

- 1.00 

- 0 . 1i7 

0 

0 

0 

0 

0 

0 

0 

9 

- 1.00 

- 0 . 1*8 

0 

0 

0 

0 

0 

0 

0 

10 

- 1.00 

- 0 . li9 

0 

0 

0 

0 

0 

0 

0 

11 

- 1.00 

- 0.50 

0 

0 

0 

0 

0 

0 

0 

12 

- 1.00 

- o . 5o 

0 

0 

0 

0 

0 

0 

0 

13 

- 1.00 

- 0.50 

0 

0 

0 

0 

0 

0 

0 

1 k 

- 1.00 

- 0.50 

0 

0 

0 

0 

0 

0 

0 

15 

- 1.00 

- 0.50 

- C.-12 

0 

0 

0 

0 

0 

0 

16 

- 1.00 

- 0.53 

- 0.19 

0 

0 

0 

0 

0 

0 

17 

- 1.00 

- 0.56 

- 0.23 

0 

0 

0 

0 

0 

0 

18 

1 

H 

• 

o 

o 

- 0.59 

- 0.25 

0 

0 

0 

0 

0 

0 

19 

- 1.00 

- 0.61 

- 0.27 

0 

0 

0 

0 

0 

0 

20 

- 1.00 

- 0.62 

- 0.29 

0 

0 

0 

0 

0 

0 

21 

- 1.00 

- 0.63 

- 0.30 

0 

0 

0 

0 

0 

0 

22 

- 1.00 

-0.6U 

- 0.31 

0 

0 

0 

0 

0 

0 

23 

- 1.00 

- 0.65 

- 0.31 

0 

0 

0 

0 

0 

0 

2b 

- 1.00 

- 0.65 

- 0.32 

0 

0 

0 

0 

0 

0 






TABLE B-VII 












k 

VALUES 





8 

9 

10 
11 
12 
13 
1 k 

15 

16 

17 

18 

19 

20 
21 
22 
23 

2k 


TABLE B-VIII 

x 1.0 06 s 1.0 


n VALUES 





-1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

-1.00 

0.50 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

-1.00 

0.25 

0.88 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

-1.00 

0.09 

0.76 

1.97 

1.00 

1.00 

1.00 

1.00 

1.00 

-1.00 

0.09 

0.61i 

0.95 

0.99 

1.00 

1.00 

1.00 

1.00 

-1.00 

0.09 

0.56 

0.93 

0.98 

1.00 

1.00 

1.00 

1.00 

-1.00 

0.09 

0.51 

0.92 

0.97 

1.00 

1.00 

1.00 

1.00 

-1.00 

0.09 

0.1*9 

0.91 

0.97 

0.99 

1.00 

1.00 

1.00 

-1.00 

0.09 

0.U7 

0.90 

0.96 

0.99 

1.00 

1.00 

1.00 

-1.00 

0.09 

0.1*6 

0.90 

0.95 

0.98 

1.00 

1.00 

1.00 

-1.00 

0.09 

0.1*6 

0.90 

0.9U 

0.98 

0.99 

1.00 

1.00 

-1.00 

0.09 

0.>45 

0.39 

0.9 h 

0.97 

0.99 

1.00 

1.00 

-1.00 

0.09 

0.1*5 

0.89 

0.93 

0.97 

0.99 

1.00 

1.00 

-1.00 

-O.lU 

0.1*5 

0.89 

0.93 

0.97 

0.99 

0.99 

1.00 

-1.00 

-0.21 

0.1*1 

0.89 

0.93 

0.96 

0.98 

0.99 

1.00 

-1.00 

-0.25 

0.37 

0.88 

0.92 

0.96 

0.98 

0.99 

1.00 

-1.00 

-0.28 

o.3lt 

0.87 

0.92 

0.96 

0.98 

0.99 

1.00 

-1.00 

-0.30 

0.32 

0.86 

0.91 

0.95 

0.98 

0.99 

1.00 

-1.00 

-0.32 

0.30 

0.86 

0.91 

0.95 

0.97 

0.99 

1.00 

-1.00 

-0.3U 

0.28 

0.85 

0.90 

0.95 

0.97 

0.99 

1.00 

-1.00 

-0.35 

0.27 

0.8U 

0.90 

0.91; 

0.97 

0.98 

1.00 

-1.00 

-0.36 

0.26 

0.8U 

0.90 

0.91* 

0.97 

0.98 

1.00 

-1.00 

-0.36 

0.25 

0.81* 

0.39 

0.91* 

0.96 

0.98 

1.00 

-1.00 

-0.37 

0.21* 

0.83 

0.89 

0.93 

0.96 

0.98 

1.00 

-1.00 

09.37 

0.21* 

0.83 

0.89 

0.93 

0.96 

0.98 

1.00 



TABLE B-VIII 









TABLE B-IX 

* i # o 06 s 2.0 


n VALUBS 


VALUES 

0 

1 

2 

3 

k 

5 

6 

7 

8 

0 

-1.00 

2.00 

2.00 

2.00 

2.00 

2.00 

2.00 

2.00 

2.00 

1 

-1.00 

1.25 

2.00 

2.00 

2.00 

2.00 

2.00 

2.00 

2.00 

2 

-1.00 

o.ee 

1.81 

2.00 

2.00 

2.00 

2.00 

2.00 

2.00 

3 

-1.00 

0.6U 

1.62 

1.95 

2.00 

2.00 

2.00 

2.00 

2.00 

k 

-1.00 

0.1*0 

1.U6 

1.92 

1.98 

2.00 

2.00 

2.00 

2.00 

5 

-1.00 

0.35 

1.33 

1.89 

1.97 

1.99 

2.00 

2.00 

2.00 

6 

-1.00 

0.26 

1.22 

1.87 

1.96 

1.99 

1.99 

2.00 

2.00 

7 

-1.00 

0.19 

l.ii* 

1.85 

1.91* 

1.98 

1.99 

1.99 

2.00 

8 

-1.00 

0.13 

1.08 

1.81* 

1.93 

1.97 

1.99 

1.99 

2.00 

9 

-1.00 

0.09 

1.03 

1.82 

1.92 

1.97 

1.99 

1.99 

2.00 

10 

-1.00 

0.05 

1.00 

1.81 

1.91 

1.96 

1.98 

1.99 

2.00 

11 

-1.00 

0.03 

0.97 

1.80 

1.90 

1.95 

1.98 

1.99 

2.00 

12 

-1.00 

0.03 

0.91* 

1.79 

1.89 

1.95 

1.98 

1.99 

2.00 

13 

-1.00 

0.03 

0.92 

1.79 

1.88 

1.91* 

1.97 

1.99 

2.00 

14 

-1.00 

0.03 

0.91 

1.78 

1.87 

1.93 

1.97 

1.98 

2.00 

15 

-1.00 

0.03 

0.90 

1.77 

1.86 

1.93 

1.96 

1.98 

2.00 

16 

-1.00 

0.03 

0.89 

1.77 

1.86 

1.92 

1.96 

1.98 

2.00 

17 

-1.00 

0.03 

0.89 

1.76 

1.85 

1.91 

1.95 

1.98 

2.00 

18 

-1.00 

0.03 

0.89 

1.76 

1.8U 

1.91 

1.95 

1.97 

2.00 

19 

-1.00 

0.03 

0.88 

1.75 

1.81* 

1.90 

1.91* 

1.97 

2.00 

20 

-1.00 

0.03 

0.88 

1.75 

1.83 

1.90 

1.91* 

1.97 

2.00 

21 

-1.00 

0.03 

0.88 

1.71* 

1.83 

1.89 

1.91* 

1.97 

2.00 

22 

-1.00 

0.03 

0.88 

1.71* 

1.82 

1.89 

1.93 

1.96 

2.00 

23 

-1.00 

0.03 

0.87 

1.71* 

1.82 

1.88 

1.93 

1.96 

2.00 

* 

-1.00 

0.03 

0.87 

1.73 

1.81 

1.88 

1.92 

1.96 

2.00 
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TABLE B-X 

^ * 1.0 OC! 3.0 


k 

VALUES 

n VALUBS 

0 

1 

2 

3 

4 

5 

6 

7 

8 

0 

-1.00 

3.00 

3.00 

3.00 

3.00 

3.00 

3.00 

3.00 

3.00 

1 

-1.00 

2.00 

3.00 

3.00 

3.00 

3.00 

3.00 

3.00 

3.00 

2 

-1.00 

1.50 

2.75 

3.00 

3.00 

3.00 

3.00 

3.00 

3.00 

3 

-1.00 

1.18 

2.50 

2.93 

3.00 

3.00 

3.00 

3.00 

3.00 

4 

-1.00 

0.97 

2.28 

2.90 

2.98 

3.00 

3.00 

3.00 

3.00 

5 

-1.00 

0.80 

2.10 

2.86 

2.96 

2.99 

3.00 

3.00 

3.00 

6 

-1.00 

0.68 

1.97 

2.83 

2.94 

2.99 

2.99 

3.00 

3.00 

7 

-1.00 

0.^8 

1.86 

2.81 

2.93 

2.98 

2.99 

2.99 

3.00 

8 

-1.00 

0.51 

1.76 

2.79 

2.91 

2.97 

2.99 

2.99 

3.00 

9 

-1.00 

o.4S 

1.71 

2.77 

2.89 

2.96 

2.99 

2.99 

3.00 

10 

-1.00 

o.4o 

1.66 

2.75 

2.88 

2.95 

2.98 

2.99 

3.00 

11 

-1.00 

0.37 

1.62 

2.74 

2.86 

2.94 

2.98 

2.99 

3.00 

12 

-1.00 

0.34 

1.58 

2.73 

2.85 

2.93 

2.97 

2.99 

3.00 

13 

-1.00 

0.32 

1.56 

2.72 

2.51* 

2.92 

2.96 

2.98 

3.00 

14 

-1.00 

0.30 

1.54 

2.71 

2.83 

2.91 

2.96 

2.98 

3.00 

15 

-1.00 

0.28 

1.52 

2.70 

2.82 

2.90 

2.95 

2.98 

3.00 

16 

-1.00 

0.27 

1.50 

2.67 

2.87 

2.89 

2.98 

2.97 

3.00 

17 

-1.00 

0.26 

1.49 

2.68 

2.80 

2.89 

2.94 

2.97 

3.00 

18 

-1.00 

0.26 

1.48 

2.67 

2.79 

2.88 

2.93 

2.97 

3.00 

19 

-1.00 

0.25 

1.47 

2.67 

2.78 

2.87 

2.93 

2.96 

3.00 

20 

-1.00 

0.24 

1.46 

2.66 

2.78 

2.86 

2.92 

2.96 

3.00 

21 

-1.00 

0.24 

1.1*6 

2.65 

2.77 

2.86 

2.92 

2.96 

3.00 

22 

-1.C0 

0.23 

1.45 

2.65 

2.76 

2.85 

2.91 

2.95 

3.00 

23 

-1.00 

0.23 

1.1*1* 

2.61* 

2.75 

2.84 

2.91 

2.95 

3.00 

24 

-1.00 

0.23 

1.1*1* 

2.61* 

2.75 

2.81* 

2.90 

2.94 

3.00 


TABLE B-X 









TABLE B-XI 

/u ■ 1.0 “ s 1*.0 



TAH.E B-XI 









TABLE B-XII 

yu » 1.0 = 5.00 


k 

VALUES 

n VALUES 

0 

1 

2 

3 

4 

5 

6 

7 

8 

0 

-1.00 

5.oo 

5.oo 

5.oo 

5.oo 

5.oo 

5.00 

5.oo 

5.oo 

1 

-1.00 

3.50 

5.oo 

5.oo 

5.oo 

5.oo 

5.oo 

5.oo 

5.oo 

2 

-1.00 

2.75 

U.62 

5.oo 

5.oo 

5.oo 

5.oo 

5.oo 

5.oo 

3 

-1.00 

2.28 

U.25 

4.90 

5.oo 

5.oo 

5.oo 

5.oo 

5.oo 

4 

-1.00 

1.95 

3.92 

a. 85 

4.97 

5.oo 

5.oo 

5.oo 

5.oo 

5 

-1.00 

1.70 

3.66 

4,79 

4.95 

4.99 

5.oo 

5.oo 

5.oo 

6 

-1.00 

1.51 

3.1*5 

4.75 

4.92 

a. 98 

4.99 

5.oo 

5.oo 

7 

-1.00 

1.37 

3.29 

4.71 

4.39 

4.97 

4.99 

4.99 

5.oo 

8 

-1.00 

1.26 

3.17 

a.68 

4.87 

4.95 

4.99 

4.99 

5.oo 

9 

-1.00 

1.17 

3.07 

a.65 

a. sa 

a. 9a 

4.98 

4.99 

5.oo 

10 

-1.00 

1.10 

2.99 

a . 63 

a .82 

4.93 

4.97 

4.99 

5.oo 

11 

-1.00 

1.05 

2.93 

a.6i 

4.60 

4.91 

4.97 

4.99 

5.oo 

12 

-1.00 

1.00 

2.88 

4.59 

4.78 

4.90 

a.96 

4.98 

5.oo 

13 

-1.00 

0.97 

2.81* 

a.58 

4.76 

a.88 

4.95 

4.98 

5.oo 

14 

-1.00 

0.91* 

2.81 

a. 56 

4.75 

4.87 

4.94 

4.97 

5.oo 

15 

-1.00 

0.93 

2.78 

a. 55 

4.73 

a.86 

4.93 

4.97 

5.oo 

16 

-1.00 

0.91 

2.76 

4.53 

4.72 

a. sa 

4.92 

4.96 

5.oo 

17 

-1.00 

0.89 

2.71* 

4.52 

4.70 

4.83 

4.91 

a .96 

5.oo 

18 

-1.00 

0.88 

2.72 

4.51 

4.69 

4.82 

4.80 

4.95 

5.oo 

19 

-1.00 

0.87 

2.71 

4.50 

a.68 

4.81 

4.89 

4.95 

5.00 

20 

-1.00 

0.86 

2.70 

4.49 

4.67 

4.80 

4.89 

4.94 

5.00 

21 

-1.00 

0.86 

2.69 

a.as 

4.65 

4.79 

4.88 

4.94 

5.00 

22 

-1.00 

0.85 

2.68 

a.a7 

a. 6a 

4.78 

4.87 

4.93 

5.oo 

23 

-1.00 

0.85 

2.67 

a.a6 

a. 63 

4.77 

a.86 

4.92 

5.00 

24 

-1.00 

0.81* 

2.66 

a.a6 

4.62 

4.76 

4.85 

4.92 

5.oo 


TABLE B-XII 









k 

n VALUBS 

VALUES 

0 

1 

2 

w 

3 

4 

5 

6 

7 

0 









1 

0.26 

0.10 

0 

0 

0 

0 

0 

0 

2 

0.75 

0.23 

o.o4 

0 

0 

0 

0 

0 

3 

0.96 

0.36 

0.10 

0.02 

0 

0 

0 

0 

k 

1.15 

0.49 

0.17 

o.o4 

0.01 

0 

0 

0 

5 

1.30 

0.63 

0.25 

0.08 

0.02 

0 

0 

0 

6 

i.4o 

0.75 

0.33 

0.12 

o.o4 

0.01 

0 

0 

7 

1.48 

0.84 

o.4l 

0.17 

0.06 

0.02 

0 

0 

8 

l.5o 

0.93 

o.49 

0.22 

0.09 

0.03 

0 

0 

9 

1.48 

0.98 

0.56 

0.28 

0.12 

o.ok 

0.01 

0 

10 

1.40 

1.01 

0.61 

0.33 

0.15 

0.06 

0.02 

0 

11 

1.30 

1.01 

0.66 

0.37 

0.19 

0.08 

0.03 

0 

12 

1.15 

0.98 

0.68 

o.4l 

0.22 

0.10 

o.oli 

0 

13 

0.96 

0.93 

0.69 

O.hh 

0.25 

0.12 

0.05 

0 

14 

0.75 

0.85 

0.69 

0.47 

0.28 

o.i4 

0.06 

0 

15 

o.5i 

0.75 

0.66 

o.48 

0.30 

0.16 

0.07 

0 

16 

0.26 

0.62 

0.62 

o.48 

0.32 

0.18 

0.08 

0 

17 

0 

0.47 

0.56 

o.47 

0.33 

0.19 

0.09 

0 

18 

-0.26 

0.32 

0.49 

0.45 

0.33 

0.20 

0.10 

0 

19 

-o.5i 

0.16 

0.41 

0.42 

0.33 

0.21 

0.10 

0 

20 

-0.75 

0 

0.31 

0.38 

0.32 

0.21 

0.10 

0 

21 

-0.96 

0 

0.21 

0.33 

0.30 

0.21 

0.11 

0 

22 

-1.15 

0 

0.17 

0.27 

0.28 

0.21 

0.11 

0 

23 

-1.30 

0 

o.i4 

0.23 

0.25 

0.19 

0.10 

0 

24 

-i.4o 


-0.34 

0.12 

0.20 
L— — 

0.22 

0.18 

0.10 

0 
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k 

VALUES 




n VALDES 




0 

1 

2 

3 

1+ 

— 

5 _ 

6 

7 

25 


-0.58 

0 

0.18 

0.20 

0.16 

0.09 

0 

26 

-1.50 

-0.70 

0 

0.11 

0.18 

0.15 

0.08 

0 

27 

-1.1(8 

-0.7l* 

0 

0.09 

O.ll* 

0.13 

0.08 

0 

28 

-1.1(0 

-0.7I* 

0 

0.07 

0.12 

0.11 

0.97 

0 

29 

-1.30 

-0.71 

0 

0.06 

0.10 

0.10 

0.06 

0 

30 

-1.15 

-0.66 

-0.21 

0.05 

0.08 

0.08 

0.05 

0 

31 

-0.96 

-0.68 

-0.29 

0 

0.07 

0.07 

o.ol* 

0 

32 

-0.75 

-0.63 

-0.33 

0 

o.ol* 

0.06 

0.01* 
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1.00 

0 

0.19 

0.1*1 

0.6U 

0.82 

0.93 

1.00 

• 

108 

• 

1 

! 

1-1*3 

0 

0.20 

0.1*2 

0.62 

0.79 

0.92 

1.00 
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APPENDIX C 


NUMERICAL SOLUTION OF A THAW- CONSOLIDATION PROBLEM 

It is intended to show by the following example a general com- 
putational procedure for solving a thaw-consolidation problem numerically# 
The physical nature of the problem and its mathematical formulation are 
discussed in Part IV* For the sake of clarity in solution technique, con- 
ditions for this example have been simplified. They may not represent a 
realistic case* 

C-01, CONDITIONS AND ASSUMPTIONS. - 

Assume for the purpose of this example that the following con- 
ditions and assumptions hold: 

(l.) The soil profile applying to this problem is the same as 
that shown in Figure 11. The sketch given below shows the profile at a 
time when the clayey silt just begins to thaw. This situation is denoted 
as zero time# 



(2.) The clayey silt has the following properties in thje un- 
frozen consolidated state: 



= 

120 

lb. per cu. ft* 

n 

s 

0.9 

at x ■ 0 

k 

- 

3.5 

—8 

x 10- ft. per seo. 

a 

V 

« 

0.1 

sq. ft. per ton 

c 

V 

• 

2.1 

-5 

x 10 sq. ft. per sec. 


It is assumed that these properties, with the exception of © do not 
change appreciably after ice segregation and during subsequent consoli- 
dation* 

(3.) Ground surface has heaved a total of 3*6", 1.6" of which 
can be accounted for by volume expansion as the soil moisture freezes# 

The remaining 2 W results from lense growth in the clayey silt which is 
assumed to be distributed uniformily with depth in the frozen zone. 

(Iu) The average rate of thawing is linear at 1/16" per hour 
and the soil follows a path similar to line BC f A of Figure 10 during 
consolidation* A discussion of this simplified case is presented in 
Part IV. 

C-02. SOIL CONDITIONS PRIOR TO THAWING. - 

The intergranular pressure p^ at the top of the clayey silt 
is approximately 0.055 tons per sq# ft. if typical values of the effective 
unit weight of gravel are used. The intergranular pressure at any depth x 
is then given by the expression: 

p T “ p + A p - 0.055 + l& . rJ&k — (i) 

* 1 2000 

- 0.055 ♦ 0.029* ..••(c- 1 ) 


C-2 


Prom condition (2.), it follows that the void ratio at any depth 
x before freezing is: 

e . = 6 , + a Ap 

Ax 1 v 

- 0.900 + 0.0029X • . . . ( C-2) 

This void ratio at any depth is indicated by point A in Figure 10. 

From condition (3*)> it may be shown that the void ratio in- 
crease Ae as a result of 2- inch, ice lens9 growth is: 

Ae= — -lllll = 0.112 ....(C-3) 

5b 

The void ratio within the frozen portion of the clayey silt immediately 
after it thaws is therefore: 

®Bx =^ 9 + e - 1.012 + C. 0029* ( C-U) 

which is equivalent to point B in Figure 10. The void ratio correspond- 
ing to point C* in the same figure is equal to 0.9095 for all depths* 

For the numerical solution, the clayey silt will be considered 
to consist of 8 layers each 8 inches thick, a subdivision identical to 
that made in Figure 11 for the hydraulic model. According to the condi- 
tion that the soil is represented by point B in Figure 10 immediately 
after thawing, the initial hydrostatic excess water pressure u is numeri- 

cally equal to p given by Equation (C-l), The excess pressure remains 

x 

equal to this value until the void ratio at any point decreases to 0,9055 
at which time the gradual transfer of overburden load from the pore water 
to the grain structure commences* 

The following table applies then to the conditions of this 


problem: 


Point 

n 

Depth 

inches 

X 

ft. 

initial 

u 

tons per 
sq. ft. 

6 x 

before 

freezing 

S x 

immediately 
after thawing 

1 

11+ 

0 

0 ( 0 . 055 ) 

0.9000 

1.0120 

2 

22 

0.67 

0 . 07 I* 

0.9019 

1.0139 

3 

30 

1.33 

0 . 091 * 

0.9039 

1.0159 

1* 

38 

2.0 

0.113 

0.9058 

1.0178 

5 

1*6 

2.67 

0.132 

0.9077 

1.0197 

6 

51+ 

3.33 

0 

0.9097 

0.9097 

7 

62 

l*.o 

0 

0.9116 

0.9116 

8 

70 

1*.67 

0 

0.9135 

0.9135 

9 

78 

5.33 

0 

0 . 9151 + 

0 . 9151 * 


C-03* NUMERICAL SOLUTION. - 

It has been demonstrated in Section 2-01-d that the finite dif- 
ference approximation of the diffusion equation for heat flow may be 
written as Equation (2 - 17c) for the one dimensional case. The corres- 
ponding equation for consolidation becomes: 

u ■ n lu _ -f- u 


n,k+l 


P (u n-i, k + Vi. k ) + (1 ' 2 ^ )u ».k •—(c-5) 


If we select a value of equal to l/3 for this solution then 
Equation (C-5) becomes: 

u = ~ (u + u + x ) ... .(C-6) 

n,k+l 3 n-l,k n,k n+l,k 


> 


Since: 




c vAt 

As 2 


....(c-7) 


and since we have selected A s = 8" then 


9 






t 



i 









At = 


I ( T2 ) 

2.1 x 1(T 5 


= 7100 sec. 

Say “2 hrs. (time interval between steps) 

Equation (C-6) applies while the soil is consolidating along a 
path given by line C*A of Figure 10. Along line BC f the excess pressure 
remains constant while the void ratio decreases. 

The finite difference approximation for void ratio change may be 


written : 


n,k+l 


n,k 


^ 'Vi, 


+ U 

.n-l,k 


^n.k^ •• ..(C-8) 


which is the comparable to Equation (2-18) of Section 2-01-f -(3 ) . 


Equation (C-8) may be rearranged 


Ka n,kn = K V + P (u n-l,k - %k> - 

and when ^ = l/3 it may be written % 

K A© = j (A^-A^) 


6 (u - U , , ) 
^ n,k n+ljk 

(C^9) 


...(C-iO) 


We must now find K before procadmg with the numerical solution. 

In each two-hour period the volume of water A.Q flowing from an 
8 -inch layer is given by Darcy ‘'s Law* 


“ Q in ~ Q out 


AH 2 AH 

* -k 7200 + k __I 7200 

= 720ac -A.H 2 ) 

= 7200k i 

As S2TTT (A“i -A u 2 ) 

* 6 X io“° (A -Au^) CU. ft. 


c-5 


L 


since i 


A® 


Av 


V _ 


A.Q 


v T(l-n) 

s 


6 x 1CT 6 (Aui - Aup) 


4 

1.9 


-6 


17.3 X 10“ (A^ -A u 2 ) 


then it follows from Equation (C-10) that* 


1 

K = T = 9*6 tons per sq. ft* 

17.3 x i0~ 6 x 2000 

The numerical solution will he carried out in a tabular form 
similar to that shown in Table II. Notation for this example is as follows 



Equations (C-6) and (C-9) govern the numerical process which is now presented. 


C-6 


•• 



KUHBUCAL SO LOTI OK OF A 
THAW - COKSOLIDATIO* PROELW 


1 

2 bn. 
1/8 in. 


8.6kOO 


> 

li bn 

lA 1 b. 


0.07k 9.733k o.07k 9.733k 


8.6kOO 


9.733k 



0.09k 9.7526 


0.U3 9.7709 


9.7526 I 0.09k I 9.7526 


9.7709 




0.132 9.7891 0.132 9.7991 



















REMARKS io flo ’ ,, occur * “BUI tha dapth of thawli* X 
panatrataa 8 Inc haw r 


63 

126 hn. 
7-7/8 1 b*. 


6k 

128 hn. 
8 Inn. 


65 

130 b 
8 - 1/8 1 


66 

132 b 

8-lA l 


105 

210 bn. 
1 3^1/8 in. 


106 

212 hn. 
13-1 A laa. 


107 

21k bn. 
1V3/8 1». 


108 

ni bn. 
1 3-1/2 la*. 


8.6kOO 0 I 8.6100 0 | 8. 6k 00 

-0.0 A -0.02k 7 -0.07k -O.03k7 




8.6kOO 


•FUJI It 


9.733k 


9«733li l 0,074 1 9 *7007 0«97L 9.681*0 0»07li | 8.7007 0.07h 8.6P20 9.021*7 8.6601 0,0002 0.6599 



Flow eonencM 
fpo« Point 2 
toward turfice 





X% la appamt 
tnat tha axcaaa 
pmaanra at 
Point 2 will to 
aaaantially dtli- 
aipatad whan tto 
thawing oana« 






127 

25k hn. 
15-7/8 In*. 


8.6kOO 


128 

256 hn. 
16 inn. 


8.6kOO 


129 

258 hn. 
16-1/8 1m. 


_2 «•«*<» 0 | 8.6kOO 

-0.0313 -0.010k -O.OklS. -0.0139 


130 

260 hn. 

14-1A 1». 


8.6582 


0 8.6582 0.0313 8.6582 0.0kl8 8.6687 

-0.09k -0.0313 -0.0627 -0.0209 -0.0522 -0.017k 


0.09k 9.7526 I 0.O9k 9.7526 0.09k 1 9.7213 0.09k 9.700k 


Flow coimikm 
fro* Point 3 
toward aurfaca 





131 

262 bn. 

16-3 A Ian. 


_0 |8.6kOO 0 | 8. 6k 00 

,Ok53 -0.0151 -O.Oh6k -0.0155 


132 

26k hn, 
16-1/2 1m, 


O .Ok53 1 8.6722 | 0.0k6k 
.0k87 -0.0162 


0.09k 9.6830 0.09k 9.6668 


Kota that th, aoll at Point 2 
■walla aa watar ^ron Point 3 
bag in* to flow onward. 


133 

266 hn. 
16-5/8 1m. 


8.6kOO 


8.6837 























































































